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THE  DYNAMICS  OF  THE  LOOP  CURRENT  AND  SHED  EDDIES  IN  A  NUMERICAL  MODEL  OF  THE 
GULF  OF  MEXICO 


HARLEY  E.  HURL6URT  and  J.  DANA  THOMPSON 

Environmental  Simulation  Branch  (Code  322),  Naval  Ocean  Research  and  Development 
Activity,  NSTL  Station,  MS  39529  USA 

ABSTRACT 

The  dynamics  of  the  circulation  in  the  Gulf  of  Mexico  have  been  investi¬ 
gated  using  simple,  efficient  numerical  models  capable  of  simulating  consis¬ 
tently  observed  dynamical  features,  including  the  Loop  Current  and  the  shedding 
of  large  anti cyclonic  eddies  from  the  Loop...  Over  150  model  experiments  were 
integrated  to  statistical  equilibrium,  typically  3-5  years. 

One  popular  hypothesis  holds  that  the  Loop  Current  sheds  anticyclonic  eddies 
in  response  to  annual  variations  in  the  inflow  through  the  Yucatan  Straits. 
However,  a  striking  result  from  the  models  is  their  ability  to  simulate  the 
observed  quasi-annual  eddy  shedding  period  with  n£  time  variations  in  the  inflow. 
The  model-predicted  eddy  diameters,  amplitudes,  and  westward  propagation  speeds 
are  also  realistic.  The  dominant  instability  mechanism  in  the  eddy  shedding  is 
a  horizontal  shear  instability  of  the  first  internal  mode,  a  barotropic  rather 
than  a  baroclinic  instability.  Therefore,  a  reduced-gravity  model  with  one 
vertical  mode  is  able  to  simulate  the  basic  dynamics  of  the  Loop  Current-eddy 
system.  Rossby-wave  theory  and  a  conservation  of  absolute  vorticity  trajectory 
analysis  were  used  to  explain  the  behavior  of  the  Loop  Current,  including  its 
northward  penetration  into  the  Gulf,  the  latitude  of  westward  bending,  the 
shedding  period  for  the  eddies,  as  well  as  their  diameter,  and  their  westward 
propagation  speed. 

A  regime  diagram  for  the  reduced-gravity  model  was  constructed  in  terms  of 
the  Reynolds  number  Re  and  the  beta  Rossby  number  Rg  =  vc/SLp  ,  where  vc  is  the 
velocity  at  the  core  of  the  current,  Lp  is  half  the  port  separation  distance  and 
3  is  differential  rotation.  Eddy  shedding  can  be  prevented  by  reducing  Re  or 
by  increasing  R3. 

Bottom  relief  acts  to  inhibit  baroclinic  instability,  yielding  solutions 
more  closely  resembling  those  from  the  reduced-gravity  model  than  the  two-layer 
flat-oottom  model.  Topography  also  influences  the  paths  of  the  shed  eddies  and, 
in  the  presence  of  sufficient  deep  water  inflow  through  the  Yucatan  Straits, 
prevents  Loop  Current  penetration,  westward  bending,  and  eddy  shedding.  In 
effect,  the  West  Florida  Shelf  acts  to  reduce  the  port  separation,  increase 
Rg,  and  shift  the  Loop  Current  into  a  stable  regime. 
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The  signatures  of  barotropic  and  barociinic  instabilities  in  the  two-layer 
Gulf  of  Mexico  model  were  studied  using  upper  and  lower  layer  pressure  fields 
and  eddy-mean  energetics.  Both  instability  processes  tend  to  drive  a  deep  flow 
characterized  by  modon  generation  and  they  exhibit  similar  vertical  phase 
relationships.  However,  in  these  experiments  the  westward  propagation  speeds 
associated  with  barociinic  instability  are  typically  two  to  three  times  faster. 


1.  I  INTRODUCTION 

Semi-enclosed  seas  are  attractive  domains  for  ocean  modeling  partly  because 
they  allow  the  use  of  numerical  grids  that  resolve  strong  meandering  currents 
and  associated  eddying  phenomena  also  found  in  major  ocean  basins.  The  Gulf  of 
Mexico  is  particularly  attractive  because  it  contains  a  major  current  system 
that  sheds  energetic  anticyclonic  eddies  which  are  comparable  in  size  to  warm- 
core  Gulf  Stream  rings.  This  system  is  illustrated  in  Fig.  1  by  the  depth  of 
the  22°C  isotherm  in  the  eastern  Gulf  of  Mexico  based  on  a  hydrographic  survey 
by  Liepper  (1970).  It  shows  the  Loop  Current  entering  from  the  south  through 
the  Yucatan  Straits  and  exiting  to  the  east  through  the  Florida  Straits.  The 
mean  transport  through  the  straits  is  -  30  nP/s  (Nowlin,  1972).  In  Fig.  1  a 
large  anticyclonic  eddy  is  about  to  break-off  from  the  Loop  Current,  as 
confirmed  by  subsequent  observations  (El liott,  1979) .  The  Loop  Current  Dene- 
trates  into  the  Gulf  and  sheds  these  large  anticyclonic  eddies  with  a  quasi¬ 
annual  period.  The  eddies  have  a  typical  radius  of  180  km  and  translate  into 
the  western  Gulf  at  a  mean  speed  of  2.4  cm/s  (Elliott,  1979). 

In  this  paper  we  present  some  basic  dynamical  ideas  and  numerical  results 
concerning  the  behavior  of  the  Loop  Current-eddy  system.  The  dynamical 
topics  include  1)  the  nature  of  the  instability  associated  with  the  eddy 
shedding,  2)  the  external  and/or  internal  factors  which  determine  the  eddy 
shedding  period,  3)  the  trajectory  dynamics  of  the  Loop  Current  and  how  they 
affect  the  penetration  of  the  Loop  into  the  Gulf,  the  eddy  shedding,  and  the 
diameter  of  the  eddies,  4)  the  existence  of  different  regimes  for  the  Loop 
Current,  5)  two  important  roles  of  topography  in  the  dynamics,  and  6)  the 
distinctive  signatures  of  barotropic  and  barociinic  instability  in  the  flow 
and  in  the  energetics.  Over  150  numerical  experiments  have  been  performed  to 
explore  the  model  parameter  space,  but  more  importantly  to  aid  in  the  formu¬ 
lation  and  testing  of  dynamical  hypotheses.  This  paper  is  both  a  distillation 
and  an  extension  of  Hurlburt  and  Thompson  (1980),  hereafter  referred  to  as  HT. 
That  paper  discusses  the  only  previous  numerical  model  of  the  Gulf  of  Mexico 
which  was  integrated  to  statistical  equilibrium  or  which  simulated  the  basic 
repetitive  features  of  the  eddy  shedding  by  the  Loop  Current. 
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Fig.  1.  Topography  of  the  22°C  isothermal  surface,  4-18  August  1966  (Alaminos 
cruise  66-A-ll)  from  Leipper  (1970).  Subsequent  data  (Elliot,  1979)  indicates 
that  an  anticyclonic  eddy  separated  from  the  Loop  Current  within  several  months 


2.  DESIGN  OF  THE  NUMERICAL  MODELS  AND  NUMERICAL  EXPERIMENTS 

The  three  numerical  models  of  HT  are  used  to  elucidate  the  Loop  Current- 
eddy  shedding  dynamics.  These  models  were  designed  to  be  as  simple  as  possible 
while  retaining  the  ability  to  simulate  the  basic  phenomena  of  interest.  The 
three  models  are  1)  reduced  gravity,  2)  barotropic,  and  3)  two-layer.  The 
first  two  are  mathematically  identical  except  for  parameter  values,  particularly 
the  gravitational  acceleration.  The  reduced  gravity  model  is  designed  to 
represent  the  first  internal  mode  and  contains  an  upper  active  layer  and  a 
lower  layer  which  is  infinitely  deep  and  at  rest.  In  the  reduced  gravity 
and  two-layer  models  the  pycnocline  is  represented  by  an  immiscible  interface 
between  two  layers  with  a  prescribed  density  contrast.  In  the  reduced  gravity 
model  bottom  topography  and  baroclinic  instability  are  not  permitted.  Two 
active  layers  is  the  minimum  to  allow  baroclinic  instability  and  to  allow 
coexistence  of  topography  and  the  pycnocline.  The  barotropic  and  reduced 
gravity  models  demonstrate  the  behavior  of  the  individual  modes  and  provide 
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insight  into  how  they  interact  in  the  two-layer  model.  They  also  allow  the 
investigation  of  some  phenomena  in  the  simplest  context. 

The  models  are  primitive  equation  on  a  B-plane  and  retain  a  free  surface. 
Using  the  hydrostatic  and  Boussinesq  approximations  and  a  right-handed 
Cartesian  coordinate  system,  the  vertically  integrated  model  equations  are 


{V  •  \vi 


•  V)  wi  +  k  x  f  W,-  =  -hiVPi 


+  Ui  -  ff i+1  )/P  +  AV2  wi 
3hi 

g~  +  V  •  Wi  =  0 

where  i  =  1,2  for  the  two-layer  model,  i=l  for  the  barotropic  and  reduced 
gravity  models  and 


(1) 

(2) 


7  '  !x  '  +  fy  ^ 


g'  =  g(p2-Pi)/p 


Pi  =  gi! 

P2  =  Pi  -  g' (h1-H1) 

A  A 

Wi  =  hiw1-  =  hi  (u1  i  +  Vi  j ) 


f  =  f0  +  3(y-y0)  (3) 

*i  =  Ti  ’  +  ^  J 


See  Appendix  A  for  symbol  definitions.  In  the  reduced  gravity  model  the  lower 
layer  momentum  equation  is  gVn  =  g'Vhj. 

Fig.  2  shows  the  model  domain  superimposed  on  a  topographic  map  of  the 
Gulf  of  Mexico.  The  20°  counter-clockwise  rotation  of  the  model  domain  is 
neglected.  The  numerical  models  were  driven  from  rest  by  prescribed  inflow 
through  the  Yucatan  Straits  (southern  port)  compensated  by  outflow  through 
the  Florida  Straits  (eastern  port).  Except  at  the  ports  the  boundaries  are 
rigid  and  in  almost  all  cases  the  no-slip  condition  is  used.  \V-j  is  prescribed 
at  the  southern  (inflow)  port.  In  most  cases  a  parabolic  inflow  profile  is 
used  forV.j.  Due  to  the  geostrophic  tilt  in  the  interface  across  the  port, 
the  velocity  maximum  is  west  of  the  center  of  the  inflow  port.  At  the  eastern 
(outflow)  port  the  normal  flow  is  self-determined  using  the  full  x-momentum 
equation.  At  inflow  points  the  boundary  condition  is  ux=0.  At  outflow 
points  the  computational  boundary  condition  is  uxx=0.  The  latter  results  in 
upstream  differencing  for  the  (Uu)x  term,  which  is  lagged  in  +ime  in  this  case. 
The  normal  pressure  gradient  is  assumed  uniform  across  the  port  and  is  determined 
by  an  Integral  constraint  requiring  the  net  outflow  through  the  eastern  port  to 
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compensate  the  inflow  through  the  southern  port.  The  tangential  velocity  com¬ 
ponent  at  the  ports  is  usually  set  at  zero  1/2  grid  distance  outside  the  physical 
domain.  This  weak  overspecification  eliminates  the  possibility  of  outflow  at 
unrealistic  angles.  Outflow  through  a  channel  modeling  the  Florida  Straits  is 
a  more  realistic  approach.  This  was  done  in  a  few  cases,  but  with  negligible 
effect.  The  semi-implicit  numerical  models  of  HT  are  used  in  this  study.  The 
implicit  treatment  of  the  external  and  internal  gravity  waves  allows  much  longer 
time  steps  in  the  numerical  integration  than  comparable  explicit  primitive  equa¬ 
tion  models  with  a  free  surface.  See  HT  for  further  discussion  of  the  numerical 
models. 


Fig.  2.  Bathymetry  of  the  Gulf  of  Mexico  based  on  U.  S.  Coast  and  Geodetic 
Survey  Chart  1007  and  soundings  on  file  at  the  Dept,  of  Oceanography,  Texas  ASM 
University.  From  Nowlin  (1972).  The  rectangle  shows  the  approximate  domain  of 
the  numerical  model.  Inflow  and  outflow  ports  are  also  indicated. 

Table  1  presents  the  parameters  of  the  pivotal  experiment  for  each  numerical 
model.  These  parameters  imply  a  maximum  upper  layer  inflow  velocity  of  70-75 
cm/s  and  an  internal  radius  of  deformation,  X  =  (g'hj)  '  /f  =  45  km,  about  four 
times  less  than  the  observed  radius  of  major  eddies  shed  by  the  Loop  Current. 

In  the  two-layer  model  the  value  of  g1  in  the  table  is  multiplied 
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by  ( H j+H£ ) /H2  to  yield  the  same  internal  values  for  the  gravity  wave  speed  as 
in  a  reduced  gravity  model.  The  inflow  transport  is  spun  up  with  a  time 
constant  of  30  days  to  minimize  the  excitation  of  high  frequency  waves.  Poten¬ 
tially  important  wind  driving  is  neglected  to  allow  focus  on  the  Loop-driven 
circulation. 

TABLE  1 


Model 

parameters  for  standard  case 

A 

lO^cn/sec"* 

8  2xl0"^cm"*sec~* 

f0 

5  x  10~5sec"1 

1  -3 

p  1  gm  cm 

g 

980  cm  sec 

*1  0 

g' 

3  cm  sec 

Ax  20  km* 

Hi 

200  m 

Ay  18. 75  km* 

h2 

2800  m 

At  1 . 5  hr 

Domain 

Size,  x.  by  y. 

1600  x  900  km 

Southern  Port  Width,  Lpw 

160  km 

Eastern  Port  Width,  Le 

150  km 

Center  of  southern  port  at  xp 

1200  km 

Center  of  eastern  port  at  y 

75  km 

r 

Upper  Layer  Inflow  Transport** 

20  x  106m3sec-1  (20  Sv) 

Lower  Layer  Inflow  Transport 

10  x  loVsec"1  (10  Sv) 

Angle 

of  inflow  from  x-axis,  Qj 

90° 

Inflow 

spin-up  time  constant 

30  days 

For  the  barotropi c  model  the  initial  maximum  depth  is  H=3000  m  and  the  inflow 
transport  is  30  Sv. 

*  for  a  given  variable 

**  also  for  the  standard  reduced  g-avity  model 


Horizontal  friction  provides  the  only  dissipation  in  the  models. 

Because  Laplacian  friction  is  a  crude  parameterization,  for  convenience 
Ah^V^v-j  was  replaced  by  AV^V^  (with  minimal  effect).  The  standard  eddy 
viscosity  (A)  is  greater  than  required  for  stable  integration  of  the  models. 
HT  showed  this  value  yields  a  constant  eddy-shedding  period  for  the  Loop 
Current.  Smaller  values  introduced  some  irregularity  into  the  period  without 
substantially  altering  the  long-term  mean.  Although  lower  eddy  viscosities 
are  utilized  in  some  experiments,  most  employ  the  larger  value  to  reduce  the 
length  of  the  integration  required  to  obtain  stable  statistics,  and  to  facili¬ 
tate  the  analysis  of  the  results.  Fig.  3  shows  the  idealized  Gulf  of  Mexico 


topography  used  in  some  of  the  numerical  experiments.  Typically  the  models 
were  integrated  five  years  to  statistical  equilibrium. 


Fig.  3.  Bathymetry  of  the  idealized  Gulf  of  Mexico  model.  The  deepest  water 
is  at  3000  m  and  the  shallowest  topography  is  400  m  deep.  The  contour  interval 
is  250  m.  (From  HT) . 

3.  AN  ATTEMPT  TO  SIMULATE  THE  EDDY  SHEDDING  BY  THE  LOOP  CURRENT 

Our  first  goal  was  to  determine  which,  if  any,  of  the  models  could  demon¬ 
strate  eddy  shedding  with  a  realistic  eddy  diameter,  amplitude,  shedding 
period,  and  propagation.  Within  the  framework  of  the  two-layer  model,  the 
first  simulation  was  made  as  realistic  as  possible  including  the  idealized 
topography  shown  in  Fig.  3.  A  longstanding  hypothesis  (Cochrane,  1965) 
maintains  that  the  Loop  Current  exhibits  an  annual  eddy  shedding  cycle  due  to 
seasonal  variations  in  the  flow  through  the  Yucatan  Straits  which  affect  the 
penetration  distance  of  the  Loop  Current.  When  the  Loop  retreats,  eddy  shedding 
is  presumed  to  occur.  Despite  this  hypothesis,  the  model  was  first  driven  by 
a  steady  inflow  to  see  if  the  Loop  Current  would  shed  eddies  due  to  purely 
internal  mechanisms.  This  might  then  establish  a  natural  frequency  for  the 
eddy  shedding. 

The  first  experiment  utilizes  the  parameters  of  Table  1  and  the  topography 
of  Fig.  3  except  that  the  upper  layer  inflow  transport  is  25  Sv,  the  lower  layer 
5  Sv.  Fig.  4  illustrates  an  eddy  shedding  cycle  from  this  experiment  using  a 
sequence  of  four  synoptic  maps  of  the  pycnocline  anomaly  (PA).  The  PA  is  the 


ower  layer  5  Sv.  (From 


deviation  of  the  interface  between  the  layers  from  its  initial  flat  elevation 
and  is  positive  downward  (upper  layer  thickness  greater  than  initial).  Fig. 

4a  shows  the  Loop  Current  has  penetrated  into  the  Gulf  and  is  beginning  to 
form  an  anticyclonic  eddy.  In  Fig.  4b  the  Loop  Current  has  bent  westward  and 
an  eddy  is  about  to  break  off.  Fig.  4c  shows  the  Loop  Current  and  an  eddy 
just  after  an  eddy-shedding  event.  In  Fig.  4d  the  eddy  has  drifted  westward 
while  the  Loop  Current  has  penetrated  further  into  the  Gulf.  Fig.  4a,  b  shows 
that  when  an  eddy  reaches  the  western  boundary,  it  drifts  northward  with 
final  decay  in  the  northwest  corner  of  the  basin. 

The  cycle  of  Loop  Current  penetration  into  the  Gulf,  westward  bending  and 
eddy  shedding  is  repeated  with  a  period  of  about  290  days,  close  to  the  quasi¬ 
annual  period  observed.  Contrary  to  the  popular  hypothesis,  time  variations 
of  the  inflow  are  not  required  for  the  model  Loop  Current  to  exhibit  realistic 
quasi-annual  eddy  shedding,  a  striking  result  first  noted  by  HT.  The  model 
also  predicts  a  realistic  eddy  diameter,  amplitude,  and  westward  propagation 
speed,  in  subsequent  experiments  HT  found  that  realistic  time  variations  in 
the  upper  layer  inflow  can  have  a  significant  influence  on  the  eddy  shedding. 
However,  the  eddy  shedding  is  dominated  by  the  natural  period,  not  the  period 
of  the  forcing,  a  topic  we  shall  not  pursue  here.  In  Section  7  we  do  examine 
an  important  effect  of  topography  and  lower  layer  inflow  through  the  Yucatan 
Straits  on  the  eddy  shedding. 

4.  A  SIMPLE  TEST  FOR  THE  INSTABILITY  MECHANISM 

The  remaining  sect'ons  are  designed  to  provide  some  insight  into  the 
dynamics  of  the  Loop  Cur,ent-eddy  system.  We  might  be  tempted  to  undertake  a 
stability  analysis  to  find  instability  mechanisms,  unstable  wavelengths,  and 
growth  rates,  but  we  would  anticipate  that  the  configuration  of  the  current 
would  be  troublesome.  However,  there  are  more  fruitful  approaches  than  this. 
Nevertheless,  we  will  start  with  one  simple  test  for  the  primary  instability 
mechanism  by  using  the  reduced  gravity  model.  If  it  produces  results  similar 
to  the  two- layer  model,  then  we  have  eliminated  baroclinic  instability  as 
an  essential  element  of  the  dynamics,  and  the  primary  instability  mechanism 
is  a  horizontal  shear  instability  of  the  internal  mode,  a  barotroDic  instability. 


1 

I 


Fig.  5  compares  (a)  the  experiment  with  topography  discussed  in  Section  3, 

(b)  a  two-layer  flat-bottom  experiment  using  the  standard  parameters  from  Table 
1,  and  (c)  a  reduced  gravity  experiment  using  appropriate  parameters  from  Table 
1.  Shown  is  a  latitude  vs.  time  (y  vs.  t)  plot  of  PA  at  a  longitude  190  km 
west  of  the  center  of  the  southern  port.  In  all  three  cases  the  PA  shows  a 
regular  progression  of  discrete  eddies  with  similar  eddy  diameter,  amplitude, 
and  shedding  period.  However,  the  experiment  with  topography  (Fig.  5a)  did  not 
begin  to  shed  eddies  for  almost  three  years,  a  point  addressed  in  Section  7. 

With  standard  parameters  the  barotropic  model  evolved  to  a  steady  state  without 
shedding  eddies,  a  matter  discussed  in  Section  6. 

The  results  shown  in  Fig.  5  lead  us  to  conclude  that  a  horizontal  shear 
instability  of  the  internal  mode  is  dominant.  The  two -layer  model  with  Fig.  3 
topography,  the  two- layer  flat-bottom  model,  and  the  reduced  gravity  model  do 
not  agree  in  all  the  parameter  space  we  explored  (see  Section  8),  but  they  do 
agree  for  a  regime  in  accord  with  observed  features  of  the  Loop  Current-eddy 
system.  Since  the  reduced  gravity  model  is  the  simplest  of  the  models  to  provide 
a  realistic  simulation,  it  is  used  in  much  of  our  analysis.  The  question  of 
barotropic  vs.  baroclinic  instability  is  addressed  further  in  Section  8  using 
eddy-mean  energetics  and  other  signatures  of  the  instability  mechanisms. 

5.  CAV  TRAJECTORIES  AND  ROSSBY  WAVE  THEORY  ELUCIDATE  THE  LOOP  CURRENT  -  EDDY 
SHEDDING  DYNAMICS 

Constant  absolute  vorticity  (CAV)  trajectories  and  Rossby  wave  theory  are 
useful  aids  in  understanding  the  dynamics  of  the  Loop  Current  and  the  eddy  shed¬ 
ding  including  the  penetration  of  the  Loop  Current  into  the  Gulf  and  the  eddy 
diameter,  shedding  period,  and  westward  propagation. 

5 . 1  CAV  trajectory  analysis 

CAV  trajectories  are  based  on  conservation  of  potential  vorticity  on  a 
8-plane  and  on  steady,  frictionless,  geostrophical ly  balanced  flow.  In  the 
reduced  gravity  model  this  implies  that  contours  of  upper  layer  depth  are 
streamlines,  and  thus  absolute  vorticity  is  also  conserved.  The  CAV 


Fiq.  5.  Time  variations  of  PA  190  km  west  of  the  center  of  the  inflow  port  for 
three  cases:  (a)  case  shown  in  Fig.  4  which  includes  topography,  (b)  standard 
two-layer  flat-bottom  case  using  parameters  from  Table  1,  and  (c)  standard 
reduced  gravity  case  using  appropriate  parameters  from  Table  1.  A  regular 
progression  of  eddies  through  the  north-south  cross  section  is  shown  in  each 
case.  The  contour  interval  is  30  m.  (From  HT). 
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r+xhc 


trajectories  are  calculated  from 


si 


co 

v  r 
c  o 


(4) 


where  0  is  the  angle  of  the  current  with  respect  to  the  positive  x-axis,  vc  is 
the  velocity  at  the  core  of  the  current,  rQis  the  radius  of  curvature,  and  the 
subscript,  0,  indicates  a  value  at  the  origin  of  the  trajectory  calculation 
(e.g.,  see  HT ;  Reid,  1972;  Haltiner  and  Martin,  1957). 

Fig.  6  shows  CAV  trajectories  superimposed  on  the  model  domain  for  vc  = 

75  cm/s,  r0=°°  and  six  different  values  of  0O.  Since  r0=«=,  the  origin  is  the 
first  inflection  point  after  inflow  (see  Fig.  4),  and  not  the  inflow  port  as 
implied  by  the  figure.  As  0Q  increases,  the  trajectory  increasingly  tends  to 
loop  back  on  itself.  When  0o=13O°,  the  CAV  trajectory  intersects  itself  at  the 
origin,  a  physically  impossible  situation  for  a  steady  flow.  Thus,  when  0O 
becomes  large  some  physical  instability  of  the  Loop  Current  can  be  anticipated. 

In  Fig.  6  00  is  varied  to  simulate  the  formation  of  an  eddy,  but  it  really 
represents  a  sequence  of  steady  state  solutions  to  (4).  Although  the  eddy- 
shedding  Loop  Current  is  not  steady,  its  evolution  is  sufficiently  slow  to 
consider  it  in  isostatic  adjustment  with  respect  to  CAV  trajectories.  In  the 
time  a  fluid  particle  in  the  Loop  moves  from  the  west  side  to  the  east  side, 
the  Loop  bends  westward  only  5  to  10%  of  the  Loop  diameter. 


5.2  Influence  of  Rossby  waves 

How  does  the  model  Loop  Current  bend  westward  when  the  angle  of  inflow 
through  the  Yucatan  Straits  is  not  varied?  We  can  gain  some  insight  into  this 
by  examining  the  continuity  equation,  (2).  If  the  mass  divergence  is  geostro- 
phic,  it  will  propagate  westward  as  a  nondispersive  internal  Rossby  wave.  (Note 
the  converse  is  not  true  near  the  equator).  For  geostrophic  divergence  (2)  becomes 


ah.  ah.  ah, 

at~ "  BVg/f  =  at~  '  cir  W  =  0 


(5) 


where  Vg  is  the  geostrophic  meridional  transport  and 

cir  =  Sg’hj/f2  =  BA2  (6) 

is  the  nondispersive  internal  Rossby  wave  speed. 

The  importance  of  nondispersive  Rossby  wave  propagation  can  be  anticipated 
from  appropriate  isolated  vortex  theory  (McWilliams  and  Flierl,  1979),  because 
1)  r/A=4  where  r  is  the  eddy  radius  and  A  is  the  internal  radius  of  deformation, 
and  2)  the  beta  Rossby  number,  RB=vc/(Br^)=!l  for  the  eddies.  From  the  linear 
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trajectori 


phase  speed  for  Rossby  waves 


cr=@/(k2+f2+X"2)  (7) 

where  k  and  l  are  zonal  and  meridional  wavenumbers,  respectively,  the  dispersive 
and  nondispersi ve  contributions  are  equal  for  circular  eddies  when  r/X=n/  V2m 
2.22.  Hence,  we  also  expect  a  significant  but  secondary  contribution  from 
dispersive  Rossby  wave  propagation. 

Fig.  7  compares  the  instantaneous  mass  divergence  (V*Wj)  and  the  geostro- 
phic  mass  divergence  (-c-j r  ^-)  for  the  standard  reduced  gravity  experiment  at 
two  different  stages  in  the  eddy  shedding  cycle.  The  usefulness  of  the  instan¬ 
taneous  mass  divergence  fields  must  be  questioned,  since  they  are  easily  domi¬ 
nated  by  rapid  oscillations  and  computational  noise.  The  absence  of  these  con¬ 
taminations  was  verified  by  showing  that  the  instantaneous  mass  divergence  and 
the  20-day  mean  are  virtually  identical.  From  (2)  the  20-day  mean  can  be 
determined  from  the  change  in  h^  in  20  days. 

Near  the  ports  and  near  the  western  boundary  the  mass  divergence  is  far 
from  geostrophic.  Note  particularly  in  Fig.  7e  the  ageostrophic  mass  conver¬ 
gence  in  the  southeastern  part  of  the  basin  associated  with  the  northward  pene¬ 
tration  of  the  Loop  Current.  Also  note  that  mass  convergence  occurs  at  the 
center  of  the  eddy  during  its  formation  {an  anticyclonic  inflow),  while  mass 
divergence  occurs  at  the  center  (anticyclonic  outflow)  after  the  eddy  separates 
from  the  Loop  and  slowly  decays.  However,  in  the  westward  bending  Loop  and  in 
the  recently  shed  eddy  the  total  mass  divergence  and  the  geostrophic  mass 
divergence  are  quite  similar.  This  clearly  demonstrates  an  important  contribu¬ 
tion  of  non-dispersive  internal  Rossby  wave  propagation  to  the  westward  bending 
of  the  Loop  Current  and  the  westward  propagation  of  the  eddies.  Other  contri¬ 
butions  such  as  nonlinear  and  dispersive  Rossby  wave  propagation  are  not 
accounted  for  here.  Still,  we  have  identified  an  important  mechanism  which 
acts  to  bend  the  Loop  Current  westward.  Thus,  it  also  produces  a  counter¬ 
clockwise  rotation  of  the  current  at  the  first  inflection  point  after  inflow 
(see  Fig.  4).  This  in  turn  produces  changes  in  the  CAV  trajectory  (Fig.  6) 
which  leads  us  to  anticipate  the  formation  and  shedding  of  an  eddy. 

5.3  Two  time  scales  associated  with  the  eddy  shedding 

Note  that  westward  bending  of  the  current  and  the  tendency  for  it  to  loop 
back  on  itself  can  be  understood  without  invoking  an  instability  mechanism. 

An  instability  mechanism  appears  essential  only  to  explain  the  separation  of 
the  eddy  from  the  Loop.  Thus,  two  time  scales  are  associated  with  the  eddy 
shedding  period:  1)  the  long  time  scale  for  the  Loop  Current  to  penetrate  into 
the  Gulf  and  bend  westward  into  an  unstable  configuration,  and  2)  the  much 
shorter  time  scale  for  the  growth  of  the  instability  as  the  eddy  separates  from 
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Fig.  7.  Shows  (a,d)  the  PA  for  the  standard  reduced  gravity  experiment, 
which  uses  the  pertinent  parameters  from  Table  1,  and  compares  the  associated 
(b,e)  instantaneous  mass  divergence  (V  •  and  (c,f)  geostrophic  mass  divergence 
(-ci r3h/3x).  Two  different  stages  of  an  eddy  cycle  are  illustrated  at  days  990 
(left)  and  1080  (right).  The  contour  intervals  are  20  m  for  PA  and  10"4  cm/s 
for  mass  divergence. 
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the  Loop  Current.  This  suggests  a  repeated  spin-up  of  the  Loop  Current  which 
eventually  becomes  unstable.  Although  it  may,  it  is  not  clear  that  the  Loop 
Current  must  satisfy  any  criterion  for  instability  during  much  of  the  eddy  cycle. 

5.4  Hypothesis  testing 

The  CAV  trajectory  analysis,  Rossby  wave  theory,  and  the  vorticity  equation 
can  be  used  to  formulate  a  number  of  quantitative  hypotheses  concerning  the 
northward  penetration  of  the  Loop  Current  into  the  Gulf,  the  latitude  at  which 
the  westward  bending  occurs,  and  the  eddy  diameter,  shedding  period,  and  west¬ 
ward  propagation.  We  will  formulate  the  appropriate  scales  and  then  test  them 
as  hypotheses  for  the  dynamics  governing  the  Loop  Current-eddy  system  predicted 
by  the  reduced  gravity  numerical  model.  The  vorticity  equation  for  the  reduced 
gravity  model  is 

3£i  3u,  3v,  „ 

iT'V  *  (ft5l>  <5T-  *  ay-1  *  6vl  ’  #7%  <8) 

where  C=v  -u  is  the  relative  vorticity. 

X  y 

5.5  Eddy  diameter.  Loop  Current  penetration,  and  latitude  of  westward  bending. 
The  beta  Rossby  number,  Rg,  is  the  ratio  of  relative  to  planetary  vorticity 

advection,  and  Rg=l  provides  a  minimum  inertial  length  scale,  Lgj=( vc/B)^2, 
over  which  6  is  important.  We  hypothesize  that  this  determines  the  latitude  at 
which  the  Loop  Current  bends  westward.  The  frictional  length  scale  over  which 
3  is  important,  Lgp=(A/8)^,  is  much  smaller.  Lgj=191  km  and  Lgp=37  km  for 
our  standard  reduced  gravity  experiment.  This  implies  that  inertia  will  prevent 
Rossby  wave  action  from  bending  the  Loop  Current  westward  at  a  higher  latitude 
than  friction. 

1/2 

We  also  find  that  Rg=l  and  r=(vc/6)  '  are  appropriate  values  for  the 
radius  of  the  eddies  formed  by  the  Loop  Current.  However,  without  further 
analysis  it  is  not  clear  why  the  eddies  from  the  Loop  Current  select  this  scale. 
McWilliams  and  Flierl  (1979)  have  studied  persistent  isolated  eddies  with 
Rg»l  and  note  that  typically  Rg >1  for  Gulf  Stream  rings. 

For  insight  into  the  scale  selection  by  the  Loop  Current  eddies,  we  turn 
to  the  CAV  trajectory  analysis  and  present  a  discussion  similar  to  HT.  Inte¬ 
gration  of  (4)  assuming  vc  =  constant  along  a  streamline  at  the  core  of  the 
current  yields 

cos©  =  cosOQ  +  5^-y2  -y/r0  (9) 

This  neglects  the  point  that  vc  f  constant  along  a  streamline,  if  there  are 
variations  in  radius  of  curvature.  The  north-south  diameter  of  an  eddy  between 
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speed  maxima  can  be  estimated  by  setting  0o=tt  at  the  southernmost  extent  (the 
origin)  and  6=0  at  the  northernmost  extent.  Then  for  r0=°°,  (9)  becomes 

d  =  2r  =  2(vc/B)1/2  (10) 

where  d  is  the  desired  diameter,  a  north-south  "dimension"  for  "stationary 
planetary  eddies"  noted  by  Rossby  (1940,  p.  82).  This  implies  that  Rg=vc/(Br^ )=1 
for  the  Loop  Current  eddies. 

The  northernmost  penetration  of  the  core  of  the  current,  b,  from  the  lati¬ 
tude  of  0O  can  be  estimated  from  (9)  by  setting  0=0  at  the  northernmost  extent. 
This  yields 


b  = 


r2  r4  2 

--  +  -2  +  2r  (1-cosO  ) 

-°  ro 


1/2 


(ID 


For  r  =“>,  (11)  reduces  to 

b  =  d  sin'20o  (12) 

Thus,  with  rQ=°°  the  maximum  amplitude  for  a  CAV  trajectory  occurs  when  0o-*m. 

This  is  what  (10)  really  represents,  since  in  this  case  the  CAV  trajectory  loops 
back  on  itself  northwest  of  the  origin  (see  Fig.  6).  For  0o=13O°  the  CAV  trajec¬ 
tory  loops  back  on  itself  at  the  origin  in  a  figure  8  (Fig.  6d)  and  sin^0o=.9. 

We  might  anticipate  that  the  horizontal  shear  instability  would  occur  when  the 
first  inflection  point  of  the  current  after  inflow  rotates  counterclockwise  to 
~  130°.  In  the  numerical  solutions  where  the  flow  is  not  steady,  potential 
vorticity  is  not  perfectly  conserved,  and  other  conditions  of  the  CAV  analysis 
are  not  perfectly  met,  we  find  that  eddy  separation  occurs  when  this  angle  is 
somewhat  >  130°  (see  Fig.  4).  In  this  case  a  CAV  trajectory  would  loop  back  on 
itself  northwest  of  the  first  inflection  point  after  inflow.  However,  a  col  is 
configured  such  that  eddy  separation  actually  occurs  southeast  of  this  inflection 
point.  From  the  standpoint  of  estimating  the  amplitude  of  the  CAV  trajectory 
(with  r0=<10)  and  thus  the  eddy  diameter,  the  value  of  0Q  when  the  eddy  separates 
from  the  Loop  Current  is  not  critical,  since  sin’200  varies  only  10?  from  0o=13O° 
to  ©0=ff. 

5.6  Tests  of  some  dynamical  hypotheses 

Tests  of  some  hypotheses  concerning  the  dynamical  behavior  of  the  Loop 
Current  are  summarized  in  Table  2.  Immediately  apparent  is  the  pervasive  role 

of  differential  rotation  8.  The  results  are  based  on  35  reduced  gravity  experi¬ 
ments  (34  for  Lgj  and  LnfJ)  from  Table  2  of  HT,  the  same  ones  they  used  in 
similar  hypothesis  testing. 
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TABLE  2 


Tests  of  some  dynamical  hypotheses  for  the  Loop  Current^eddy  system. 


%  bias 

Standard  reduced  gravity  case 
correlation  numerical  theoretical 

1. 

Eddy  radius 

r=(vc/S)1/2  6 

.87 

186  km 

191  km 

2. 

Distance  from  the  southern  boundary  to  the 

latitude  of  westward  bending  by 

the  Loop  Current. 

Lgl  =  ( vc/e) 1/2  =•  r  -7 

.75 

201  km 

186  km 

-1 

.77 

201  km 

191  km 

3. 

Maximum  northward  penetration  of 

the  Loop 

Current. 

Lnp  =  Lei  +  b  =  3r  -2 

4 

.99 

574  km 

560  km 

.89 

574  km 

573  km 

4. 

Westward  propagation  speed  of  the  eddies. 

cir  =  BA2  =  Bg'hj/f2  40 

.99 

3.21  cm/s 

4.57  cm/s 

cr  =  B/(k2+£2+A"2)  -2 

.97 

3.21  cm/s 

3.49  cm/s 

5. 

Eddy  shedding  period. 

1- 

pe  =  A0  +  Ajr(l+cos0j )  J/ce 

.95 

327  days 

359  days 

- 

.96 

327  days 

338  days 

The  %  bias  and  the  linear  correlation  are  statistics  for  the  theoretical  pre¬ 
diction  vs.  the  values  observed  in  the  reduced  gravity  numerical  model.  They 
are  based  on  34  to  35  numerical  experiments  from  Table  2  of  HT,  the  same  ones 
they  used  in  their  hypothesis  testing  for  similar  quantities.  The  %  bias  e 
((m0-mp)/mp)  where  m0  and  mp  are  the  means  of  the  observed  and  predicted  values, 
respectively.  The  two  rightmost  columns  present  the  results  observed  and 
predicted  for  the  standard  reduced  gravity  experiment  which  uses  the  pertinent 
parameters  from  Table  1.  Results  from  two  tests  are  presented  for  Lgj,  Lnp 
and  pe.  On  the  upper  line  values  of  r  and  ce  observed  in  the  numerical  model 
were  used  in  the  predictor.  On  the  lower  line  the  theoretical  values  were  used 
for  r  and  ce  (i.e.  cr  for  ce). 

In  estimating  r,  Lgj,  and  Lnp  the  maximum  speed  at  inflow  was  used  for  vc. 
One-half  the  north-south  diameter  between  speed  maxima  was  used  for  the  eddy 
radius  from  the  reduced  gravity  numerical  model.  The  distance  from  the  southern 
boundary  to  the  southern  end  of  the  eddy  diameter  was  used  for  Lgj  and  from  the 
southern  boundary  to  the  northern  end  of  the  diameter  for  Lnp.  These  were 
measured  as  the  eddy  center  passed  a  longitude  110  km  west  of  the  western 
boundary  of  the  inflow  port.  This  was  also  close  to  the  inflection  point  which 


18 


exhibited  a  large  angle  at  this  time  {see  Figs.  4  &  7).  Thus,  the  theoretical 
estimate  of  b  was  simplified  by  setting  00=ir  and  r0=^>.  This  yields  b=2r. 

The  agreement  between  the  theoretical  estimates  of  r,  L,,j,  and  Lp  and  the 
values  calculated  from  the  numerical  model  is  remarkably  good.  In  most  cases 
the  agreement  for  r  was  within  5  .  However,  at  low  Reynolds  numbers  and  low 
latitudes,  (10)  overestimated  the  model  value  by  >  10%.  In  the  experiments  with 
low  Reynolds  numbers  downstream  attenuation  of  the  current  appears  to  explain 
this.  At  low  latitudes  the  assumption  of  isostatic  adjustment  with  respect  to 
CAV  trajectories  is  not  as  good  due  to  an  increased  westward  propagation  speed. 
Also  r/A  is  less,  so  the  Loop  and  the  eddies  are  more  subject  to  dispersion. 

In  estimating  L^j  and  Lnp)  values  of  r  from  the  numerical  model  were  used  for 
the  upper  line  and  theoretical  values  for  the  lower  line. 

In  estimating  A,  he  from  Table  2  of  HT  was  used  for  h^  and  the  value  of  f 
at  the  latitude  of  the  eddy  center  was  used.  The  correlation  between  c-r  and 
the  westward  propagation  speed  of  shed  eddies  is  extremely  high,  but  the  theo¬ 
retical  speed  is  40  greater  than  observed  in  the  reduced  gravity  numerical 
model.  When  the  dispersive  contribution  is  included,  the  mean  Rossby  wave 

speed  and  the  mean  observed  westward  propagation  speed  differ  by  <  2%.  The 
2  2 

value  of  k  +t  was  estimated  by  assuming  circular  eddies  and  using  the  theore¬ 
tical  values  for  r. 

In  Table  2  a  regression  equation  is  used  to  test  the  hypothesis  that  the 
eddy  shedding  period  is  a  multiple  of  the  time  required  for  an  eddy  to  move  one 
eddy  radius  westward.  The  multiple  depends  on  the  angle  of  inflow  for  the 
current,  Oj.  The  upper  line  presents  the  results  when  values  measured  from  the 
numerical  model  are  used  for  r  and  cg  (the  westward  propagation  speed).  On 
the  lower  line  theoretical  values  (cr  and  c0  and  (10)  for  r)  are  used.  The 
regression  coefficients  are  A0=45.1  days  and  Aj=4.67  for  the  upper  line  and 
A0=-.2  days  and  Aj=5.34  for  the  lower  line.  The  regression  coefficient,  ]2Aj  , 
is  a  kind  of  inverse  Strouhal  number  (nondimensi onal  period)  if  we  take  the 
eddy  diameter,  d,  as  the  appropriate  diameter,  and  the  westward  propagation 
speed  of  the  eddies  (rather  than  the  injection  velocity)  as  the  appropriate 
velocity.  The  Strouhal  numbers  (S=2/Aj)  implied  by  the  two  regression  results 
are  .43  for  the  upper  line  and  .37  for  the  lower  line.  If  d  was  the  half¬ 
wavelength  for  a  continuous  wavetrain,  then  the  Strouhal  number  would  be  S=.5. 

The  eddy  shedding  period  also  exhibits  a  weak  dependence  on  the  eddy 
viscosity  A,  primarily  because  this  affects  the  amount  of  entrainment  or 
decrainment  downstream.  The  result  is  an  eddy-shedding  period  which  increases 
with  increasing  Reynolds  number,  Re=vc^L/A  where  vci  is  the  maximum  velocity 
at  inflow  and  L  is  half  the  inflow  (southern)  port  width.  At  high  Reynolds 
numbers,  secondary  circulations  of  the  Loop-eddy  system  become  significant 
and  introduce  some  irregularity  into  the  eddy  shedding  period. 
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6.  REGIMES  FOR  THE  LOOP  CURRENT  IN  THE  REDUCED  GRAVITY  MODEL 

6. 1  The  eddy- shedding  regime  (E) 

The  preceding  section  examined  some  of  the  dynamics  of  the  eddy  shedding 
by  the  Loop  Current  and  demonstrated  the  usefulness  of  CAV  trajectory  analysis 
and  Rossby  waves  in  explaining  the  behavior.  It  also  demonstrated  the  important 
role  of  differential  rotation  (8)  in  most  aspects  of  the  eddy-shedding  dynamics. 
We  will  call  the  eddy-shedding  regime  the  E  regime. 

6.2  The  steady  westward  spreading  regime  (W) 

This  section  investigates  the  existence  of  other  flow  regimes  in  the 
neighboring  parameter  space.  One  such  regime  was  found  by  reducing  the  Reynolds 
number,  Re  =  vcL/A,  where  L  is  the  half-width  of  the  port.  This  acts  both  to 
damp  physical  instabilities  in  the  current  and  to  decrease  potential  vorticity 
conservation,  thus  reducing  the  tendency  of  the  current  to  loop  back  on  itself. 
The  result  is  a  steady  solution  with  a  westward  bending  Loop  Current  as  shown 
in  Fig.  8a.  We  will  call  this  the  W  regime.  Fig.  8a  was  obtained  by  increasing 
the  eddy  viscosity  (A)  from  10^  to  3xl0^cm^/s  in  the  standard  reduced  gravity 
model.  Steady  linear  viscous  solutions  (where  8vj=Av^c^)  and  weakly  nonlinear 
solutions  belong  to  the  W  regime.  As  mentioned  in  Section  5,  the  latitude  of 
westward  bending  after  inflow  is  determined  by  the  larger  of  L^j  and  Lgp.  As 
noted  by  HT,  the  mean  of  the  standard  reduced  gravity  experiment  over  an  eddy 
cycle  is  very  similar  to  Fig.  8a.  This  implies  that  in  the  mean  the  eddies  drive 
a  northward-flowing  western  boundary  current. 

Between  the  E  and  W  regimes  there  is  a  transition  regime  (T)  with  eddies 
superimposed  on  a  westward  bending  Loop  Current.  This  is  illustrated  in  Fig.  8b 
using  the  standard  reduced  gravity  model  except  that  A=2.5  x  10  cm  /s.  Experi¬ 
ments  that  exhibited  both  eddy  shedding  and  an  unbroken,  westward  bending  zero 
contour  (like  Fig.  8b)  were  assigned  to  the  T  regime. 

6.3  Steady  source-sink  regime  (N) 

A  third  major  regime  we  call  the  N  regime  is  found  by  increasing  the  beta 
Rossby  number,  Rb=vc/(6L^),  where  Lp  is  half  the  distance  separating  the  centers 
of  the  inflow  and  outflow  ports.  The  N  regime  occurs  when  Rg>2.  At  sufficiently 
low  Reynolds  numbers  the  transition  between  the  N  and  W  regimes  is  determined 
by  the  beta  Ekman  number,  Eg=A/(BLp  ),  as  discussed  in  the  next  subsection.  In 
either  case  the  N  regime  is  characterized  by  a  steady  source-sink  flow  with  no 
westward  bending  by  the  Loop  Current.  It  is  illustrated  in  Fig.  9  by  two  vari¬ 
ations  on  the  standard  reduced  gravity  experiment.  In  Fig.  9a  B=0  and  Rg=®. 

Note  that  8>0  is  essential  for  the  westward  bending  of  the  Loop  Current  and  the 
eddy  shedding  shown  in  Figs.  4  and  7.  With  Re  sufficiently  high,  we  would  expect 
some  instability  to  occur  in  the  source-sink  current,  but  not  the  quasi-annual 
eddy  shedding  exhibited  by  the  Loop  Current. 


Fig.  a.  Illustration  of  (a)  steauy  regime  W  and  (b)  time-dependent  regime  T. 

The  standard  parameters  for  the  reduced  gravity  model  were  used  except  that  in 
(a)  A  =  3  x  107  cm^/sec  and  in  (b)  A  =  2.5  x  297  cm^/s.  The  contour  interval  is 
20  m.  In  regime  W  the  Loop  Current  bends  westward  and  is  steady.  In  regime  T 
eady  shedding  is  superimposed  on  a  westward  bending  loop  instead  of  being 
discrete  as  in  Fia.  7a,  d. 


Fig.  9.  Illustration  of  the  steady  regime  N  for  the  standard  reduced  gravity 
experiment  except  that  g  =  0  and  (b)  a  section  of  land  was  added  to  approximate 
the  west  Florida  shelf.  Steady  mixed  regime  M  is  shown  in  (c)  with  parameters 
chosen  to  approximate  those  of  Mellor  and  Blumberg  (1981),  including  Re  =  10.8 


In  a  second  experiment  the  port  separation  was  reduced  by  inserting  a  land 

mass  which  approximates  the  location  of  the  West  Florida  Shelf  (Fig.  9b).  In 

this  case  Lp  was  measured  using  a  port  centered  75  km  south  of  the  western 

boundary  of  the  inserted  land  mass.  From  Section  5  r=L  T=(v  /rO  thus  the  beta 

2  i  c 

Rossby  number  can  also  be  expressed  as  Rg=(r/Lp)  or  Rg=(L:,[/Lp)  .  Since  the  N 
regime  occurs  for  Rg>2,  this  implies  a  critical  port  separation  21  -  *\/2r  =  V2L, 
From  the  definition  of  Rg,  both  H  and  Lp  play  a  similar  role  in  determining  the 
N  regime.  Otherwise,  their  roles  are  not  similar.  As  long  as  Rg  is  small 

enough  for  eddy  shedding  to  occur,  the  E  regime  is  quite  insensitive  to  Lp,  but 

most  aspects  of  the  dynamics  are  very  sensitive  to 

The  N  regime  occurs  when  the  Loop  Current  reaches  a  steady  state  before 
penetrating  far  enough  into  the  Gulf  to  bend  westward.  From  (12)  and  Fig.  6d  a 
CAV  trajectory  with  rQ=“>  can  penetrate  into  the  Gulf  with  no  westward  bending 
a  distance  b  =  ^2r  =V2L.j  (the  same  as  the  critical  port  separation).  Using 
this  as  an  estimate  of  the  maximum  possible  Loop  Current  penetration  in  the  N 
regime  is  consistent  with  our  present  numerical  results  (see  HT  for  further 
discussion). 

The  transition  between  regimes  N  and  E  is  quite  abrupt  and  a  transition 
regime  has  not  been  found  in  any  of  our  numerical  experiments.  The  transition 
between  regimes  N  and  W  is  broader.  In  this  mixed  regime  (M)  the  westward 
spreading  of  the  Loop  Current  is  significantly  less  than  it  would  be  if  the  port 
separation  were  infinite.  The  numerical  simulation  of  the  Loop  Current  by 
Blumberg  and  Mel  lor  (1981)  belongs  to  this  M  regime.  Fig.  9c  shows  a  reduced 
gravity  analog  of  that  experiment  using  our  estimates  of  appropriate  parameters, 
including  Re=10.8  and  Rg= 1.19  (See  Appendix  B  for  other  parameters  used). 
Considering  our  neglect  of  the  20°  counterclockwise  rotation  of  the  basin,  our 
pattern  for  the  Loop  Current  is  remarkably  similar  to  that  from  the  much  more 
complicated  and  expensive  model  of  Blumberg  and  Mellor  (1981).  Experiments  were 
assigned  to  the  M  regime  if  the  central  contour  on  inflow  exhibited  any  westward 
bending,  but  the  amplitude  of  the  Loop  PA  was  <  30  of  the  maximum  a  distance 
(vcl-/6)2  west  of  the  center  of  the  inflow,  where  vc1-  is  the  maximum  velocity  at 
i nflow. 

6 . 4  Stabi lity  regime  diagram 

Three  nondimensional  parameters  from  the  vorticity  equation  (Eq.  (8))  play 
an  important  role  in  determining  the  stability  regimes  for  the  Loop  Current  in 
the  reduced  gravity  and  flat-bottom  barotropic  models.  They  are  the  Reynolds 
number  (Re),  the  beta  Rossby  number  ( Rg ) ,  and  the  beta  Ekman  number  (Eg).  Only 
two  of  these  are  independent.  Provided  the  same  scales  are  used,  Eg=Rg/Re. 

Fig.  10  shows  the  parameter  space  occupied  by  the  various  regimes  on  a  stability 
diagram  of  Re  vs.  Rg-  Eddy  shedding  occurs  for  Rg<2  and  Re>25  for  the  reduced 


1  tnree  stability  regimes  ana  transition  regimes  between  them  for  a  reduce* 
e  the  Reynolds  number  (Re)  and  the  beta  Rossby  number  (Rg).  The  regimes  , 
eady  westward  spreading  {Fig.  8a),  (N)  steady  source-sink  (Fig.  9a,  b),  ( 


gravity  model,  Re>40  for  the  flat-bottom  barotropic  model  (not  shown).  Apparently, 
the  dispersive  nature  of  Rossby  waves  in  the  barotropic  model  is  a  catalyst  which 
leads  to  a  higher  critical  Reynolds  number  (Rec).  The  transfer  of  energy  to 
the  lower  layer  in  the  two-layer  models  seems  to  have  a  similar  effect  on  Rec. 

At  low  Reynolds  number  the  transition  between  regimes  N  and  W  depends  upon 
Eg=RgL/(ReLp) .  The  correction  factor  L/Lp  is  required  because  ye  have  not  used 
Lp  in  our  definition  of  Re.  This  transition  is  clearly  defined  on  our  stability 
diagram  only  if  L/ Lp  is  fixed.  From  Section  5.5  Rg  and  Eg  can  be  expressed  in 
terms  of  the  latitude  of  westward  bending  by  the  Loop  Current,  i.e.  RR=(LoT/Ln)2 
and  Eg=(Lgp/Lp)  .  Whether  Rg  or  Eg  determines  the  transition  to  the  N  regime 
depends  upon  the  relative  importance  of  L^p,  the  frictional  length  scale,  and 
Lpj,  the  inertial  length  scale,  in  determining  the  latitude  of  westward  bending 
by  the  Loop  Current.  From  the  regime  diagram,  the  transition  between  these  two 
criteria  occurs  when  2.5  Lgp=Lgj/  V2. 

Other  parameters  not  accounted  for  here  such  as  inflow  angle,  basin 
geometry,  and  basin  orientation  may  also  have  some  influence  on  the  regime 
selection.  Additional  influences  in  the  two-active-layer  model  are  discussed 
in  the  next  two  sections. 

7.  PREVENTION  OF  EDDY  SHEDDING  BY  TOPOGRAPHY  AND  DEEP-WATER  INFLOW  THROUGH 
THE  YUCATAN  STRAITS 

In  this  section  the  two-layer  model  is  used  to  demonstrate  a  dramatic 
effect  of  topography  on  the  eddy-shedding  by  the  Loop  Current  when  there  is 
sufficient  deep-water  inflow  through  the  Yucatan  Straits.  Fig.  11a  shows  the 
domain-averaged  upper  layer  (upper  curve)  and  lower  layer  kinetic  energy  vs. 
time  for  a  two-layer  experiment  the  same  as  shown  in  Fig.  4,  except  that  during 
the  first  six  years  of  model  integration  the  lower  layer  inflow  was  zero.  This 
includes  the  topography  of  Fig.  3.  At  the  beginning  of  year  7  the  lower  layer  inflow 
was  increased  to  10  Sv  with  a  time  constant  of  30  days.  While  the  lower  layer 
inflow  was  zero,  eddy  shedding  occurred  in  a  manner  similar  to  that  shown  in 
Fig.  4.  The  signature  of  the  eddy  shedding  cycle  is  depicted  in  the  upper  layer 
energy  curve  of  Fig.  11a.  When  the  lower  layer  inflow  was  increased,  the  eddy 
shedding  ceased  and  the  solution  evolved  to  a  steady  state  as  shown  in  Fig.  11a. 

This  steady  solution  is  shown  in  Fig.  lib,  c  in  terms  of  the  PA  (Fig.  lib)  and 
the  lower  layer  pressure,  P2  (Fig.  11c).  The  PA  depicts  a  source-sink  flow  like 
regime  N,  while  the  flow  in  the  lower  layer  follows  the  f/h  contours  of  the  topo¬ 
graphy.  HT  uses  a  kinematic  analysis  to  illuminate  the  dynamics  of  this 
phenomenon.  From  the  continuity  equation  (2),  the  divergence  term  hjV  •  \vj  is 
balanced  by  the  advective  term  Wj*Vhj  in  a  steady  state.  Since  the  flow  is 
nearly  geostrophic,  \vi*7hj  -  Wjg-Vhi  =  w^g’Vhj,  where  w^g  is  the  geostrophic 
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Fig.  11.  Results  for  the  bottom  topography  experiment  shown  in  Fig.  4  except 
that  Tj  =  25  Sv  and  T2  =  0  until  year  6  when  T2  increases  to  10  Sv.  (a)  domain- 
averaged  upper  layer  (top  curve)  and  lower  layer  kinetic  energy,  (IC  =  1.5).  (b) 
Nearly  steady  PA  at  day  2880.  The  contour  interval  is  20  m.  (c)  Lower  layer 
pressure  normalized  by  density  ( P2 )  at  day  2880.  The  contour  interval  is  .25  m^/s^ 


velocity  component  in  layer  i.  The  magnitude  of  Wgg'Vhj  is  greatest  when 
relatively  strong  lower  layer  currents  flow  at  large  angles  to  contours  of  hj. 
Comparison  of  Fig.  lib  and  11c  shows  that  this  occurs  where  the  Loop  Current 
intersects  a  current  following  the  f/h  contours  of  the  West  Florida  Shelf.  If 
the  lower  layer  current  is  strong  enough  for  the  advection  to  balance  the 
divergence  in  (2)  associated  with  the  approaching  Loop,  then  the  interface  deepening 
and  Loop  Current  penetration  in  this  region  are  halted. 

Thus,  the  West  Florida  Shelf  along  the  eastern  boundary  of  the  domain  (see 
Fig.  3)  and  lower  layer  flow  act  in  conjunction  to  effectively  reduce  the  port 
separation  by  locally  limiting  the  northward  penetration  of  the  Loop  Current. 

If  this  results  in  Rg>2,  when  we  measure  2Lp  as  the  distance  between  the  center 
of  the  inflow  port  and  the  point  where  the  upper  and  lower  layer  currents  inter¬ 
sect,  then  the  upper  layer  current  exhibits  the  source-sink  flow  characteristic 
of  the  N  regime  described  in  Section  6.  This  is  illustrated  by  comparing  Fig.  lib 
with  Fig.  9b,  the  reduced  gravity  experiment  with  a  land  mass  in  the  location  of 
the  West  Florida  Shelf. 

In  the  experiment  shown  in  Fig.  11  the  Loop  Current  had  already  penetrated 
far  into  the  Gulf  and  shed  eddies  when  the  lower  layer  inflow  was  increased. 

When  this  increase  occurred,  a  current  following  the  f/h  contours  developed  in 
the  lower  layer  and  the  advection  term  (w^Vhj)  began  to  exceed  the  divergence 
term  in  the  continuity  equation  where  the  Loop  Current  crossed  the  shelf  slope, 
causing  h^  to  decrease  there.  Thus,  the  Loop  Current  retreated  southward  until 
an  equilibrium  occurred  near  the  southern  end  of  the  shelf. 

In  terms  of  vorticity  dynamics,  the  northward  penetration  of  the  Loop 
Current  is  halted  when  the  interaction  between  the  topography  of  the  West 
Florida  Shelf  and  the  pressure  field  results  in  a  near  balance  between  the 
pressure  torques  and  the  nonlinear  terms  in  the  mass  transport  vorticity 
equation.  HT  discuss  and  document  this  topic  in  more  detail. 

The  results  in  this  section  suggest  that  certain  time  variations  in  the 
deep  flow  through  the  Yucatan  Straits  may  have  a  greater  effect  on  the  Loop 
Current  than  fluctuations  in  the  upper  ocean  current  through  the  Strait.  The 
effects  of  the  latter  are  discussed  by  HT  but  notin  this  paper. 
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8.  BAROTROPIC  VS.  BAROCLINIC  INSTABILITY  AND  THE  IMPORTANT  ROLE  OF  TOPOGRAPHY 

Section  6  examined  stability  regimes  for  the  reduced  gravity  model.  In 
Section  7  we  studied  a  steady  regime  for  the  two-layer  model  which  results  from 
a  particular  combination  of  topography  and  deep  flow  following  f/h  contours. 

In  this  section  we  will  investigate  additional  regimes  of  the  two  active  layer 
model,  confining  ourselves  to  unsteady  regimes  with  eddies.  In  particular  we 
will  search  for  cases  of  barotropic,  baroclinic,  and  mixed  instability. 

The  reduced  gravity  and  barotropic  models  can  exhibit  barotropic,  but  not 
baroclinic  instability.  The  reduced  gravity  model  has  demonstrated  an  eddy 
shedding  regime  with  a  horizontal  shear  instability  of  the  first  internal  mode 
(a  "barotropic"  instability)  which  produces  a  remarkable  simulation  of  observed 
features  of  the  Loop  Current  -  eddy  shedding  system.  We  know  from  Section  4 
that  in  some  cases  the  two  active  layer  model  exhibits  similar  results  with 
similar  parameters.  In  this  section  we  will  investigate  the  potential  importance 
of  baroclinic  instability  and  the  role  of  topography  in  determining  its  importance. 

8. 1  Eddy-mean  energetics 

We  begin  this  investigation  by  surveying  the  eddy-mean  energetics  for  the 
seven  numerical  experiments  listed  in  Table  3.  We  then  illustrate  some  charac¬ 
teristic  features  of  the  different  regimes  using  synoptic  maps  of  upper  and  lower 
layer  pressure  (pi  and  P2)  and  curves  of  domain-averaged  energy  vs.  time.  Fig.  12 
shows  the  eddy-mean  energetics  in  terns  of  energy  box  diagrams.  Fig.  12a  labels 
the  energy  transfers.  See  Appendix  A  for  symbol  definitions  and  Appendix  C  for 
the  energy  transfer  integrals.  All  of  the  model  domain  was  used  in  calculating 
the  energetics  except  the  parts  within  100  km  of  the  eastern  boundary  and  37.5 
km  of  the  southern  boundary.  Thus,  the  eastern  and  southern  boundaries  of  the 
energetics  calculations  are  open.  Kinetic  energy  and  pressure  work  fluxes 
through  these  open  boundaries  are  represented  by  arrows  at  the  top  (bottom)  of 
the  Kj  (Kg)  boxes.  In  all  cases  most  of  the  energy  flows  into  K^.  In  some 
cases  there  is  significant  efflux  from  K|,  but  always  must  less  than  the  Kj-*Kj 
transfer.  The  arrows  pointing  outward  from  the  sides  represent  dissipation  of 
a  particular  type  of  energy  due  to  Laplacian  horizontal  friction.  Arrows  between 
the  boxes  represent  conversions  of  energy  from  one  type  to  another  as  indicated 
by  the  direction  of  the  arrow. 

Fig.  12b  shows  the  eddy-mean  energetics  for  Experiment  1  on  Table  3,  a  two- 
layer  flat-bottom  experiment  using  the  standard  parameters  given  in  Table  1. 

The  energy  conversion  is  characteristic  of  a  barotropic  instability.  The 

potential  energy  transfer  is  actually  reversed  with  eddy  potential  energy  (P1) 
feeding  the  mean. 

Angling  the  inflow  27°  west  of  normal  in  the  standard  two-layer  flat-bottom 
model  (Experiment  2)  produced  a  dramatic  change  in  the  eddy-mean  energetics  which 
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TABLE  3 


Model 

experiments  discussed  In  Section  8 

Exp  f1 

Differences  from  standard  two-layer 
flat-bottom  experiment  in  Table  1 

Figures  in  HT  from  these 
experiments  (HT  Fig.  #) 

1 

None 

Figs.  5b,  lib,  14a 

2 

0 j  =  27°,  Sv2  =  0 

Figs.  12,  13a,  14d 

3 

Reduced  gravity,  Yucatan  and 

Florida  Straits  added  to  model 
domain 

Figs.  20a,  21a, 

HT  Table  2  case  RG32 

4 

A  =  3  x  106  cm2/s,  Sv2  =  0 

Experiment  in  HT  Figs.  10, 
11a,  14c  exhibits  similar 
behavior 

5 

A  =  3  x  105  cm2/s ,  Sv2  =  0 

Fig.  3  topography 

6 

Svj  =  25,  Sv2  =  0 

Figs.  24b,  25b 

7 

Svj  =  25,  Sv2  =  0, 

Fig.  3  topography 

Figs.  24c,  25c 

Sv.  is 

the  inflow  in  layer  i  in  106  m3/s  or  Sv. 

is  shown  in  Fig.  12c.  (The  lower  layer  inflow  was  also  reduced  to  zero,  but 
other  experiments  show  this  has  a  relatively  minor  role  in  altering  the  energetics 
in  this  case).  Fig.  12c  illustrates  a  classic  signature  of  baroclinic  instability 
in  the  eddy-mean  energetics  with  P-+P'  dominating  the  mean  to  eddy  energy 
transfer  and  feeding  the  upper  and  lower  layers  almost  equally.  There  is  even 
a  reverse  cascade  in  the  kinetic  energy  (Kj+Kj)  with  eddies  feeding  the  mean 
flow. 

The  eddy-mean  energetics  for  the  reduced  gravity  model  (Experiment  3  and 
Fig.  12d)  illustrates  a  pure  barotropic  instability,  since  this  model  excludes 
baroclinic  instability.  In  this  case  the  dominant  mean  to  eddy  energy  transfer 
is  K->K' .  Even  though  this  is  a  pure  barotropic  instability,  there  is  a  net 
transfer  from  P-*P'.  Thus,  the  existence  of  such  a  transfer  does  not  necessarily 
imply  any  contribution  from  baroclinic  instability. 
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Fig.  Id.  Eddy-mean  energetics  for  the  experiments  listed  in  Table  3:  (a)  labels 

for  the  energy  pathways  and  energy  reservoirs,  (b)  energetics  for  Experiment  1, 
standard  flat-bottom  case,  (c)  Experiment  Z,  flat  bottom  with  non-normal  Inflow, 
(d)  Experiment  3,  reduced  gravity,  (e)  Experiment  4,  flat  bottom  with  low 


viscosity,  (f)  Experiment  5,  same  as  the  preceding  but  with  topoaraphy,  (g) 

Exoeriment  6,  flat  bottom,  (h)  Experiment  7,  same  as  preceding  but  with  topography 

15 

The  energy  reservoirs  are  in  units  of  10  joules  and  the  energy  transfers  are 
in  units  of  108  joules/s. 
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Fig.  12e  is  particularly  interesting  because  it  illustrates  a  mixed  insta¬ 
bility  and  because  it  demonstrates  the  value  of  separating  the  kinetic  energy 
into  upper  and  lower  layer  components.  These  results  were  obtained  primarily 
by  reducing  the  eddy  viscosity  in  the  two-layer  flat-bottom  model  by  a  factor  of 
three  (Experiment  4).  If  Kj  and  l<2  were  combined  to  produce  a  4-box  diagram, 
the  results  would  look  much  like  those  for  the  reduced  gravity  model  and  we 
might  conclude  that  this  is  a  case  of  barotropic  instability.  In  contrast,  the 
6-box  diagram  (Fig.  12e)  illustrates  a  striking  result.  Although  a  barotropic 
energy  conversion  (Kj-H<J)  is  dominant  in  the  upper  layer,  the  lower  layer  eddies 
are  fed  almost  equally  by  transfers  from  P-*-P‘  and  Kj-*-P‘.  In  view  of  the  reduced 
gravity  results,  this  is  insufficient  evidence  for  an  important  contribution  from 
baroclinic  instability  to  the  lower  layer  eddies.  Additional  evidence  for  this 
will  be  provided  shortly. 

Figs.  12e  and  12f  compare  the  results  for  Experiments  4  and  5.  The  experi¬ 
ments  are  identical  except  that  Experiment  4  (Fig.  12e)  has  a  flat  bottom  and 
Experiment  5  (Fig.  12f)  includes  the  idealized  Gulf  of  Mexico  topography  shown 
in  Fig.  3.  The  topography  strongly  suppresses  any  baroclinic  instability.  With 
the  topography  added,  the  energy  box  diagram  (Fig.  12f)  indicates  a  strong  baro- 
topic  instability  (Kj+Kj)  and  a  strong  reverse  potential  energy  flux  (P'+P). 

Figs.  12g  and  12h  again  compare  experiments  with  and  without  the  topography 
of  Fig.  3  (Experiments  6  and  7).  They  differ  from  the  preceding  by  a  three-fold 
increase  in  the  eddy  viscosity  and  a  25%  increase  in  the  upper  layer  inflow. 

The  experiment  with  the  topography  (Fig.  12h)  exhibits  essentially  the  same 
energy  pathways  as  the  previous  frame  with  the  same  topography  (Fig.  12f).  The 
reverse  potential  energy  transfer  (P'+P)  is  even  stronger.  Almost  1/3  of  the 
eddy  energy  makes  a  complete  circuit.  Although  this  reverse  transfer  is  clearly 
augmented  by  the  topography,  it  is  not  restricted  to  experiments  with  topography 
(see  Fig.  12b).  Without  the  benefit  of  the  energetics  analysis,  HT  correctly 
identified  the  baroclinical ly  unstable  case  (Fig.  12c)  and  a  mixed  instability 
case  similar  to  Fig.  12e.  However,  they  also  conjectured  that  Experiment  6 
(Fig.  12g)  is  a  case  of  mixed  instability.  This  is  not  corroborated  by  the 
domain  integrated  eddy-mean  energetics.  Eddies  in  both  layers  are  fed  by  energy 


conversions  appropriate  for  a  barotropic  instability,  the  lower  layer  fed 
indirectly  via  energy  transfer  from  the  upper  layer.  In  this  case,  increasing 
the  eddy  viscosity  has  suppressed  the  contribution  from  baroclinic  instability. 
Experiment  6  (Fig.  12g)  differs  from  Experiment  1  (Fig.  12b)  by  having  25% 
greater  inflow  in  the  upper  layer  and  none  in  the  lower  layer.  The  energy 
pathways  in  the  two  experiments  are  similar,  but  Fig.  12g  shows  more  energy 
transfer  to  the  lower  layer  and  lacks  the  reverse  (P'-*P)  transfer  of  the  standard 
flat-bottom  experiment  (Fig.  12b). 

The  energy  transfers  in  all  these  experiments  are  strongly  inhomogeneous  in 
space.  Thus,  as  stressed  by  Harrison  and  Robinson  (1978),  energy  transfers 
averaged  over  the  model  domain  may  not  be  characteristic  of  any  important 
subregion. 


8. 2  Kinetic  energy  vs.  time 

In  the  following  discussion,  we  will  illustrate  features  of  the  flow  which 
are  characteristic  of  the  three  regimes  identified  in  the  eddy-mean  energetics 
with  barotropic,  baroclinic,  and  mixed  instabilities.  We  will  utilize  the  four 
experiments  which  do  this  most  simply  and  clearly,  (a)  Experiment  1  for  barotropic 
instability  with  a  flat  bottom,  (b)  Experiment  7  for  barotropic  instability  with 
topography,  (c)  Experiment  2  for  baroclinic  instability,  and  (d)  Experiment  4 
for  mixed  instability.  Fig.  13  shows  the  curves  of  K1  and  K2  vs.  time  for  these 
four  experiments. 

Fig.  13a,  b  represents  the  barotropically  unstable  experiments  and  clearly 
shows  a  relatively  long  period  for  the  eddy  shedding  cycle,  273  days  for  Experi¬ 
ment  1  (Fig.  13a)  and  250  days  for  Experiment  7  (Fig.  13b).  Fig.  13c  shows  a 
much  faster  57  day  oscillation  for  the  baroclinically  unstable  experiment. 
Experiment  2.  The  corresponding  reduced  gravity  experiment  (not  shown),  in 
which  baroclinic  instability  is  not  permitted,  has  a  284  day  period.  The  period 
in  Fig.  13c  is  very  similar  to  that  found  by  Holland  and  Lin  (1975)  for  mid¬ 
latitude  mesoscale  eddies  in  a  two-layer  model  with  baroclinic  instability. 

They  also  noted  a  similar  maximum  in  Kj  near  the  onset  of  baroclinic  instability 
which  is  followed  by  a  rise  in  K2.  We  have  not  found  this  type  of  signature  in 
any  of  our  barotropically  unstable  experiments.  Fig.  13d  shows  Kj  and  K2  vs. 
time  for  Experiment  4,  the  experiment  for  which  the  eddy-mean  energetics  (Fig.  12e) 
suggest  a  mixed  instability.  Two  periods  which  are  not  harmonically  related  are 
clearly  indicated,  a  long  period  of  300  days,  which  is  typical  of  barotropically 
unstable  experiments,  and  a  much  shorter  56  day  period  similar  to  that  for  the 
baroclinically  unstable  experiment.  Also  notable  are  the  dramatic  spikes  in 
K2  lagging  the  maximum  in  Kj .  This  resembles  the  behavior  of  K2  at  the  onset 
of  baroclinic  instability  shown  in  Fig.  13c. 
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Fig.  13.  Average  kinetic  energy  over  the  rectangular  domain  (upper  curve  for 
upper  layer)  for  (a)  Experiment  1,  the  standard  two-layer  flat-bottom  case, 

(b)  Experiment  7,  with  the  topography  of  Fig.  3, (c)  Experiment  2,  with  non-normal 
inflow,  and  (d)  Experiment  4,  Identical  to  Experiment  1  but  with  A  =  3  x  10°cm  /s 
and  Svg  =  0.  The  value  of  IC  is  (a)  .5,  (b)  1.0,  (c)  .55,  (d)  .7. 
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8.3  Modon  generation  in  the  barotropical ly  unstable  experiments 

We  begin  examining  the  characteristic  features  of  the  flow  in  different 
regimes  by  studying  two  experiments  where  the  eddy-mean  energetics  indicate 
barotropic  instability.  One  experiment  has  a  flat  bottom,  and  the  other  includes 
the  idealized  Gulf  of  Mexico  topography  shown  in  Fig.  3.  The  two  experiments 
which  illustrate  the  basic  features  of  this  flow  in  the  simplest  and  clearest 
fashion  are  Experiment  1  with  Fig.  12b  energetics  (the  standard  flat-bottom 
experiment)  and  Experiment  7  with  topography  and  Fig.  12h  energetics. 

Fig.  14  shows  synoptic  views  of  pj  and  p2  for  Experiment  1.  At  day  1710, 
p^  (Fig.  14a)  shows  the  Loop  Current  penetrating  into  the  basin  and  beginning 
to  bend  westward.  An  eddy  shed  earlier  lies  in  the  western  Gulf.  A  characteristic 
feature  of  the  barotropical ly  unstable  experiments  is  the  generation  of  a  modon 
in  the  lower  layer  as  the  Loop  Current  begins  to  form  an  eddy  (Fig.  14b).  The 
relationship  between  p^  and  p£  can  be  seen  clearly  by  superimposing  the  fields. 

The  modon  intensity  tends  to  follow  that  of  the  generating  eddy  in  the  upper 
layer.  The  axis  of  the  modon  is  oriented  close  to  the  direction  of  propagation 
by  the  upper  layer  vortex  with  the  anticyclonic  member  leading  and  the  cyclonic 
member  trailing.  The  orientation  of  the  modon  generated  here  is  quite  different 
from  that  found  by  McWilliams  and  Flier!  (1979)  for  isolated,  nearly  circular 
vortices,  but  the  tendency  of  the  eddy  in  the  upper  layer  to  propagate  toward 
the  member  of  the  modon  with  like  rotation  is  similar.  However,  in  this  case 
the  westward  propagation  speed  of  the  modon  slightly  exceeds  that  of  the  upper 
layer  vortex.  Thus  the  flow  actually  becomes  more  baroclinic  and  in  Fig.  14a,  b 
we  see  the  anticyclonic  eddy  in  the  upper  layer  situated  over  the  cyclonic  member 
of  the  modon.  This  behavior  is  common  but  not  universal  in  our  numerical  experi¬ 
ments.  It  is  quite  unlike  the  coupled  behavior  of  the  isolated  baroclinic 
vortex  and  barotropic  modon  studied  by  McWilliams  and  Flierl  (1979).  In  their 
results  the  modon  member  with  rotation  unlike  the  baroclinic  vortex  eventually 
broke  away  and  the  barotropic  and  baroclinic  vortices  tended  to  become  super¬ 
imposed  and  to  approach  a  state  of  deep  compensation  (no  signature  of  the  vortex 
in  the  lower  layer).  When  the  upper  layer  vortex  reaches  the  western  boundary 
and  propagates  northward  (Fig.  14c),  it  is  again  associated  with  a  modon  in  the 
lower  layer  (Fig.  14d)  and  again  the  modon  is  oriented  in  the  direction  of 
propagation  with  the  like  (anticyclonic)  member  leading  and  the  opposite  member 
trailing. 

Fig.  14b,  d  shows  an  additional  interesting  phenomenon  which  occurs  in  the 
lower  layer.  If  the  lower  layer  were  integrated  separately  as  a  barotropic 
model,  the  solution  would  evolve  to  a  steady,  westward-bending  Loop  Current  as 
in  regime  W.  The  solution  would  be  similar  to  Fig.  8a  but  the  loop  would  bend 
westward  at  a  lower  latitude  because  { vc/6 )^  is  much  less  (see  Section  5.5). 


The  lower  layer  inflow  velocity  at  the  core  of  the  current  is  only  3.35  cm/sec. 


Experiment  1.  The  contour  interval  is  .5  m  /s  for 


In  this  experiment  with  two  active  layers,  the  lower  layer  flow  is  quite  different 
from  the  barotropic  prediction.  During  the  formation  of  the  modon,  the  Loop 
Current  in  the  lower  layer  exhibits  a  source-sink  flow  (Fig.  14b).  The  modon 
is  the  dominant  flow  and  prevents  the  natural  westward  bending  of  the  Loop. 

When  the  modon  moves  away  (Fig.  14d),  the  lower  layer  Loop  bends  far  to  the  west 
and  sheds  a  weak  eddy  almost  in  phase  with  the  upper  layer  (Fig.  14c).  In  the 
mean  (Fig.  15)  both  the  upper  and  lower  layers  exhibit  a  westward-bending  loop 
like  the  W  regime.  As  predicted  by  (vc/fs)'J  and  the  higher  inflow  velocity  in 
the  upper  layer,  the  mean  loop  in  the  upper  layer  bends  westward  at  a  higher 
latitude.  The  lower  layer  also  exhibi ts  counter-rotating  and  zonally  elongated 
mean  gyres  north  of  the  loop.  These  are  driven  by  downward  flux  of  eddy  energy 
from  the  upper  layer  and  may  also  influence  the  latitude  of  the  mean  Loop  Current 
in  the  lower  layer. 

Fig.  16  shows  a  synoptic  view  of  p^  and  P2  at  day  1760  for  a  barotropically 
unstable  experiment  with  Fig.  3  topography  (Experiment  7  with  Fig.  12h  energetics). 
This  experiment  exhibits  coupled  upper  layer  vortex,  lower  layer  modon  behavior 
similar  to  the  flat  bottom  experiment,  except  that  the  modon  is  mostly  confined 
to  the  abyssal  plain.  Another  difference  is  that  the  upper  layer  vortex  remains 
between  the  modon  pair.  The  modon  is  partially  steered  by  the  topography. 
Apparently,  the  back  interaction  from  the  modon  to  the  upper  layer  is  sufficient 
that  the  trajectory  of  the  upper  layer  vortex  is  also  modified  by  the  topography. 
Fig.  17  compares  upper  layer  eddy  trajectories  for  Experiments  6  and  7,  two 
experiments  with  no  flow  through  the  ports  in  the  lower  layer.  The  experiments 
are  identical  except  that  Experiment  6  has  a  flat  bottom  and  Experiment  7  includes 
Fig.  3  topography.  Because  Experiment  7  includes  no  flow  through  the  ports  in 
the  lower  layer,  there  is  no  current  following  the  f/h  contours,  unlike  Fig.  11c. 
The  addition  of  such  a  current  had  no  major  effect  on  the  modon,  provided  the 
current  was  weak  enough  to  permit  the  normal  eddy  shedding  to  occur. 

Although  eddy  activity  in  the  lower  layer  modified  the  propagation  of  the 
upper  layer  vortex,  the  propagation  of  both  the  upper  layer  vortex  and  the 
associated  modon  was  dominated  by  internal  Rossby  wave  propagation  in  both  the 
reduced  gravity  and  two  active  layer  experiments  which  exhibited  discrete  eddy 
shedding  and  a  horizontal  shear  (barotropic)  instability  of  the  internal  mode. 

As  we  will  see  shortly,  this  is  not  the  case  in  the  experiment  with  a  baroclinic 
instability. 
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Fig.  16.  (a)  p-i  and  (b)  Op  at  day  1760  for  Experiment  7.  The  contour  interval 

2?  ?  2 
is  .5  m  /s  for  p1  and  .05  m  /s  for 


8.4  Flow  characteristics  associated  with  baroclinic  instability 

Fig.  18  shows  a  synoptic  view  of  p^  and  p£  for  Experiment  2  where  the  upper 
layer  inflow  is  angled  27°  west  of  normal  and  there  is  no  flow  through  the  ports 
in  the  lower  layer.  The  eddy-mean  energetics  (Fig.  12c)  indicate  the  occurrence 
of  baroclinic  instability.  This  experiment  exhibits  modon-like  generation  similar 
to  that  earlier  associated  with  a  barotropic  instability.  However,  the  eddies 
tend  to  be  smaller  and  the  greater  population  of  eddies  tends  to  mask  the  modon 
character  of  the  eddy  generation.  The  upper  and  lower  layer  eddies  near  the 
eastern  part  of  the  Loop  bear  a  phase  relationship  which  is  similar  to  the 
barotropically  unstable  experiments.  The  modon  axis  is  oriented  close  to  the 
direction  of  propagation  of  the  anticyclonic  eddy  in  the  upper  layer,  with  the 
anticyclonic  modon  member  leading  and  the  cyclonic  one  trailing.  One  difference 
is  that  the  modon  axis  is  south  of  the  upper  layer  vortex.  Thus,  the  lower  layer 
eddies  tend  to  be  strongest  under  the  westward-flowing  arm  of  the  Loop  Current  as 
expected  for  a  baroclinic  instability  (Gill,  et  al  ,  1974;  Philander,  1976). 

Later,  in  the  central  basin  the  leading  modon  member  shifts  northward,  away  from 
the  westward  propagating  vortex  in  the  upper  layer.  The  trailing  vortex  remains 
under  the  westward  branch  of  the  Loop.  Thus  the  modon  axis  is  no  longer  aligned 
with  the  direction  of  propagation.  In  general,  the  lower  layer  eddies  tend  to 
be  elongated  meridionally  in  the  eastern  part  of  the  basin  where  they  originate 
and  zonal ly  in  the  western  part  of  the  basin  where  they  decay.  In  the  western 
part  of  the  basin  the  eddies  also  show  some  tendency  toward  barotropy.  Except 
for  the  initial  meridional  elongation,  these  tendencies  are  consistent  with 
results  presented  by  Rhines  (1977). 

The  most  dramatic  difference  between  the  experiments  with  barotropic  and 
baroclinic  instability  lies  in  the  propagation  speed  of  the  eddies.  In  the 
barotropically  unstable  experiments  with  discrete  eddies  the  internal  Rossby 
wave  speed  associated  with  the  upper  layer  vortex  exerts  primary  control  on  the 
propagation  in  both  layers.  Even  though  the  eddies  in  the  baroclinically 
unstable  case  are  smaller,  they  propagate  westward  at  ~  10  cm/sec,  typically 
2  to  3  times  faster  than  in  the  barotropically  unstable  experiments.  Although 
it  is  difficult  to  estimate  an  appropriate  shear  velocity,  the  propagation  speeds 
in  our  numerical  model  are  quite  consistent  with  those  for  baroclinic  instability 
in  a  linearized  two-layer  model  with  a  horizontally  uniform  basic  flow(Pedlosky,  1979). 
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Fig.  18.  (a)  pi  and  (b)  p?  at  day  2450  for  Experiment  2.  The  contour  interval 
P  o  1  L  2  2 

is  .5  m  /s  for  p^  and  .05  m  /s  for  p2. 
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Gill,  et  al  (1974)  have  suggested  the  upper  to  lower  layer  phase  shift  as 
a  means  of  detecting  baroclinic  instability.  In  our  results  we  find  this  is  not 
very  useful  because  the  barotropic  instability  which  occurs  in  the  upDer  layer 
generates  a  modon  in  the  lower  layer  with  upper-lower  layer  phase  relationships 
which  are  much  like  those  of  the  baroclinic  instability.  In  our  results  the 
westward  propagation  speed  of  the  eddies  is  a  much  clearer  distinguishing 
characteristic. 

8.5  Flow  characteristics  of  a  mixed  instability 

Fig.  19  shows  two  synoptic  views  of  p^  and  p£  for  Experiment  4  which  has  a 
flat  bottom,  no  inflow  in  the  lower  layer,  and  one-third  the  eddy  viscosity  of 
the  experiments  discussed  in  Sections  8.3  and  8.4.  The  eddy-mean  energetics 
(Fig.  12e)  suggest  that  a  mixed  instability  occurs  in  this  experiment.  Since 
there  is  no  flow  through  the  ports  in  the  lower  layer,  all  the  energy  in  the 
lower  layer  is  received  from  the  upper  layer.  Apart  from  this  the  flows  in  the 
two  layers  are  much  more  independent  than  those  discussed  in  the  two  preceding 
subsections. 

Eddies  in  the  lower  layer  propagate  westward  at  approximately  the  external 
Rossby  wave  speed  (~  10  cm/sec),  and  those  with  like  rotation  pass  a  given  point 
with  a  periodicity  of  about  60  days.  Unlike  the  experiments  discussed  in 
Sections  8.3  and  8.4  there  is  no  clear  phase  relationship  between  the  eddies 
in  the  lower  layer  and  the  eddy  which  forms  on  the  Loop  Current  in  the  upper 
layer.  This  is  true  during  most,  but  not  all,  of  the  eddy-shedding  cycle  of  the 
Loop  Current.  In  this  experiment  the  eddy-shedding  period  is  about  300  days  and 
is  depicted  as  a  slow  oscillation  in  Kj  (Fig.  13d).  There  is  a  back  interaction 
from  the  lower  layer  eddies  to  the  Loop  Current  in  the  upper  layer  which  causes 
a  strong  undulation  of  the  Loop  with  approximately  a  60  day  period.  This  is 
depicted  in  Fig.  13d  as  the  high  frequency  oscillation  in  Kp  Except  for  this 
undulation,  the  Loop  Current  penetrates  into  the  Gulf,  bends  westward  and  begins 
to  form  an  eddy  structure  just  as  In  the  barotropical ly  unstable  experiments,  but 
near  the  time  an  eddy  would  break  off  (Fig.  19a)  in  a  barotropical ly  unstable 
experiment,  something  quite  different  occurs.  The  Loop  Current  suddenly  shoots 
far  to  the  west  at  a  speed  appropriate  for  baroclinic  instability  and  breaks 
into  a  series  of  smaller  eddies.  During  this  process  lower  layer  eddies  under 
the  south  side  of  the  Loop  strengthen  dramatically  and  the  upper  and  lower  layer 
eddies  develop  distinct  phase  relationships.  An  anticyclonic  eddy  in  the  lower 
layer  leads  the  westward  advance  on  the  Loop  Current.  At  this  stage  phase 
relationships  in  the  upper  and  lower  layers  are  very  similar  to  those  for  the 
baroclinically  unstable  case,  and  they  exhibit  the  same  differences  from  the 
barotropically  unstable  experiments.  These  phase  relations  disintegrate  as  soon 
as  the  rapid  westward  advance  of  the  Loop  Current  is  halted.  Thus  we  have  a 


picture  of  episodic  baroclinic  instability  associated  with  a  small  part  of  each 
300  day  eddy-shedding  cycle  of  the  Loop  Current.  This  instability  is  strong  enough  to 
show  up  in  the  domain-averaged  eddy-mean  energetics  (Fig.  12e)  and  to  provide 
a  sharp  spike  in  the  curve  of  «2  vs.  time  (Fig.  13d).  Day  1720  (Fig.  19b)  is 

near  the  foot  of  the  last  spike  and  day  1760  (Fig.  19d)  is  near  the  top  of  it. 

The  much  weaker  coupling  for  the  two  layers  than  found  in  either  the 

barotropically  or  baroclinically  unstable  experiments  is  explained  in  part  by 

the  peculiar  episodic  nature  of  the  baroclinic  instability  in  Experiment  4  and 
in  part  by  the  3  times  lower  eddy  viscosity.  Because  the  eddies  in  the  lower 
layer  are  governed  by  external  Rossby  wave  propagation,  they  are  dispersive  in 
nature.  With  the  lower  eddy  viscosity  they  are  not  dissipated  as  soon  after 
generation  and  have  greater  opportunity  to  disperse  and  fill  the  basin.  The 
importance  of  dispersion  in  spreading  the  eddy  population  in  the  lower  layer  has 
been  noted  by  Rhines  (1977)  and  in  the  study  of  isolated  vortices  by  McWilliams 
and  Flierl  (1979). 

A  comparison  of  Fig.  16  (for  the  experiment  with  idealized  Gulf  of  Mexico 
topography)  and  Fig.  19  (for  the  flat-bottom  experiment  with  a  mixed  instability) 
indicates  how  the  topography  can  suppress  the  episodes  of  baroclinic  instability 
found  in  the  latter  case.  When  the  topography  of  Fig.  3  is  present,  the  eddies 
in  the  lower  layer  are  mostly  confined  to  the  abyssal  plain.  Lower  layer  eddy 
generation  over  the  strongly  sloping  topography  is  prevented  because  the  eddy 
flow  would  have  to  cross  the  closely  packed  f/h  contours  at  large  angles, 
behavior  not  anticipated  in  geostrophically  balanced  flow  which  conserves 
potential  vorticity.  The  strong  eddies  in  Fig.  19d  which  form  under  the  westward¬ 
flowing  branch  of  the  Loop  Current  lie  over  the  region  of . the  Campeche  Bank  and 
the  slope  of  the  bank.  Thus,  they  are  prevented  from  forming  when  the  topography 
of  Fig.  3  is  included.  We  can  now  appreciate  why  the  two-layer  model  with  topo¬ 
graphy  produces  results  more  like  the  reduced  gravity  model  than  does  the  two- 
layer  flat-bottom  model.  If  the  westward  branch  of  the  Loop  Current  were  to 
flow  over  the  abyssal  plain,  we  might  expect  the  model  to  exhibit  episodes  of 
baroclinic  instability  even  when  the  topography  is  included. 
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9.  SUMMARY  AND  CONCLUSIONS 

In  this  paper  we  have  demonstrated  the  ability  of  simple  numerical  models  to 
perform  remarkable  simulations  of  the  Loop  Current  -  eddy  system  in  the  Gulf  of 
Mexico.  Tne  models  were  designed  for  computational  efficiency  and  simplicity, 
retaining  only  the  essential  physics  and  characteristics  of  the  Gulf  required  to 
simulate  the  basic  dynamical  behavior  of  the  system.  The  simplicity  of  the  models 
facilitated  the  analysis  of  the  system  dynamics  and  the  computational  efficiency 
allowed  us  to  perform  over  150  multi-year  integrations.  The  efficiency  also 
allowed  us  to  use  horizontal  resolution  adequate  to  investigate  interesting 
regions  of  the  parameter  space,  regions  where  time  dependent  eddies  dominated 
the  circulation. 

The  numerous  numerical  experiments  explored  the  model  parameter  space  and 
aided  in  the  formulation  and  testing  of  dynamical  hypotheses.  The  most  salient 
results  from  the  model  are  summarized  in  the  seven  points  which  follow. 

1)  The  simple  models  were  able  to  simulate  the  anticyclonic  eddy  shedding  by 
the  Loop  Current  and  to  simulate  eddies  with  realistic  diameters,  amplitudes, 

and  westward  propagation  (Figs.  1,4).  Most  striking  was  the  ability  of  the  models 
to  simulate  the  observed  quasi-annual  period  of  the  eddy  shedding  with  rm  time 
variations  in  the  inflow  through  the  Yucatan  Straits  (southern  port).  This  is 
contrary  to  the  popular  hypothesis  that  the  Loop  Current  sheds  eddies  in  response 
to  annual  variations  of  the  inflow. 

2)  The  reduced  gravity  model  proved  to  be  the  simplest  model  able  to  simulate 
the  basic  dynamics  of  the  Loop  Current  and  the  eddy  shedding  (Fig.  5).  This 
indicates  that  baroclinic  instability  is  not  an  essential  element  of  the  dynamics. 
Instead  a  horizontal  shear  instability  of  the  first  internal  mode  (a  barotropic 
instability)  is  the  dominant  instability  mechanism. 

3)  We  have  demonstrated  the  usefulness  of  CAV  trajectory  analysis  and  internal 
Rossby  waves  in  explaining  the  eddy-shedding  behavior  of  the  Loop  Current,  including 
the  eddy  diameter,  the  Loop  Current  penetration  into  the  Gulf,  the  latitude  of 
westward  bending  by  the  Loop  Current,  the  westward  speed  of  eddy  propagation,  and 
the  eddy  shedding  period  (see  Table  2).  The  role  of  differential  rotation  (B) 

is  pervasive.  This  theory  also  showed  that  it  is  not  necessary  to  invoke  an 
instability  mechanism  to  explain  the  westward  bending  of  the  Loop  Current,  nor 
the  tendency  for  it  to  loop  back  on  itself.  This  suggests  that  an  instability 
may  be  essential  only  to  explain  the  final  eddy  separation  from  the  Loop  Current. 

4)  Two  steady  regimes  were  found  in  the  parameter  space  neighboring  the  eddy 
shedding  regime.  The  parameter  space  occupied  by  each  regime  was  depicted  on  a 
regime  diagram  for  the  reduced  gravity  model  (Fig.  10). 

5)  In  the  presence  of  sufficient  deep  water  inflow  through  the  Yucatan  Straits, 
the  Florida  Shelf  topography  may  prevent  Loop  Current  penetration,  westward  bending 
and  eddy  shedding  by  effectively  reducing  the  port  separation.  This  shifts  the 
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Loop  Current  into  one  of  the  stable  regimes  (Section  7).  Certain  time  varia¬ 
tions  in  the  deep  flow  through  the  Yucatan  Straits  may  have  a  greater  effect 
on  the  Loop  Current  than  fluctuations  in  the  upper  layer  flow  through  the  strait. 

6)  Bottom  topography  plays  another  important  role  by  inhibiting  baroclinic 
instability.  Thus  the  reduced  gravity  solutions  were  more  like  two-layer  solutions 
with  the  idealized  Gulf  of  Mexico  topography  (Fig.  3)  than  two-layer  flat-bottom 
solutions.  The  topography  also  demonstrated  some  ability  to  steer  eddies  in  the 
upper  ocean  through  back  interaction  from  eddies  in  the  lower  layer  which  were 
mostly  confined  to  the  abyssal  plain. 

7)  Finally,  we  examined  the  characteristic  signatures  of  barotropic  and 
baroclinic  instability  in  the  pressure  fields  of  both  layers  and  in  the  eddy- 
mean  energetics.  In  both  cases  there  was  a  tendency  for  eddies  in  the  upper 
layer  to  drive  a  modon  in  the  lower  layer.  The  upper  and  lower  layer  phase 
relations  were  surprisingly  similar  for  both  types  of  instability,  but  the  west¬ 
ward  propagation  speeds  associated  with  baroclinic  instability  were  typically 
two  to  three  times  faster. 
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APPENDIX  A 
List  of  Symbols 


Hj,  H2(x,y) 
hi,  h2 


horizontal  eddy  viscosity 
nondispersive  internal  Rossby  wave  speed 
internal  Rossby  wave  speed  including  dispersion 
beta  Ekman  number,  A/(BLp  ) 

Coriolis  parameter;  f0  taken  at  southern  boundary  (y0) 

acceleration  due  to  gravity 

reduced  gravity,  g(p2  -  p^/p 

initial  thicknesses  of  the  layers 

instantaneous  local  thickness  of  the  layers 


i 


p  _ 

p,  p,  p' 


g 

VI  i ,  w2 
x,  y,  z 

V  *l 

Ax,  Ay 
6 


P »  Pj »  P2 

Tx  Ty 

Ti '  Ti 


2  2 

kinetic  energy  Js(u  +  v  ) ;  of  the  mean  flow;  mean  of  the 
eddy  flow,  respectively,  for  layer  i. 

zonal  and  meridional  wave  numbers,  respectively 

halfwidth  of  the  southern  port 

minimum  frictional  length  scale  over  which  3  is  important, 

(A/6)1/3 

minimum  inertial  length  scale  over  which  8  is  important, 

(vc/6)% 

maximum  northward  penetration  of  the  Loop  Current 

half  the  port  separation  distance 

2  2 

potential  energy  *sp(gnj  +  9^2 )  i  the  mean  flow;  mean  of 
the  eddy  flow,  respectively. 

upper  layer  density-normalized  pressure,  gnj. 

lower  layer  density-normalized  pressure,  gnj  -  g'lhj-Hj) 

eddy  shedding  period 

2 

beta  Rossby  number,  v  /6L 

L  p 

Reynolds  number,  vcL/A 

eddy  radius 

time 

time  increment  in  the  numerical  integration 

x-directed  components  of  current  velocity 

maximum  inflow  speed 

speed  at  the  core  of  the  current 

geos  trophic  meridional  transport 

hj  Wj,  h2  w2 

tangent  plane  Cartesian  coordinates:  x  positive  eastward, 
y  positive  northward,  z  positive  upward 

east-west  and  north-south  domain  size 

horizontal  grid  increments 

differential  rotation,  df/dy 

relative  vorticity  v  -  uu 

free  surface  anomaly;  height  of  the  free  surface  above 
its  initial  uniform  elevation;  n^  =  +  h2  -  Hj  -  H2 

n2  =  Hi  +  ni  "  =  h2  -  H2  =  -PA 

angle  of  inflow  with  respect  to  the  positive  x-axis 

internal  radius  of  deformation 

densities  of  oea  water 

x  and  y  directed  tangential  stresses  at  the  top  (i)  and 
bottom  (i  +  1)  of  layer  i 
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APPENDIX  B 


New  Reduced  Gravity  Experiments  for  the  Regime  Diagram  (Fig.  10) 


Experiment 

vci 

L 

A 

LP 

6 

Re 

rb 

Regime 

155 

54.4 

80 

4.0 

227.1 

2.0 

10.9 

.527 

W 

157 

54.1 

80 

4.0 

150.9 

2.0 

10.8 

1.19 

M 

158 

73.4 

80 

2.0 

203.5 

2.0 

29.4 

.89 

E 

159 

73.2 

80 

2.5 

203.5 

2.0 

23.4 

.88 

T 

160 

40.3 

80 

7.3 

203.5 

.49 

4.4 

1.99 

N 

161 

41.4 

80 

1.6 

203.5 

.537 

20.7 

1.86 

M 

162 

41.2 

80 

1.6 

203.5 

.69 

20.6 

1.44 

W 

163 

40.5 

80 

4.6 

203.5 

1.21 

7.0 

.81 

M 

164 

40.8 

80 

4.6 

203.5 

.69 

7.1 

1.43 

N 

units:  vci 

in  cm/s 

;  L, 

Lp  in  km; 

A  in  10 

7  cm2/s 

;  8  in 

10*13  cm' 

^sec'1 

See  Apoendix  A  for  symbol  definitions.  Note  that  L  and  Lp  are  the  half-port 
width  and  half-port  separation,  respectively.  The  other  parameters  are  the 


same  as  in  Table  1  except  that  (1)  in  Experiments  155  and  157  a  land  mass  was 
inserted  in  the  location  of  the  West  Florida  Shelf  (See  Fig.  9b,  c)  and  the 
center  of  the  eastern  port  was  taken  to  be  75  km  south  of  the  western  end  of 
this  land  mass,  (2)  the  southern  port  was  centered  at  xp  =  1000  km  in  Experiment 
155  and  at  xp  =  1160  km  in  Experiment  157,  and  (3)  the  inflow  transoort  was 
14  Sv  in  Experiments  155  and  157,  and  10  Sv  in  Experiments  160-164. 

APPENDIX  C 

Derivation  of  Eddy-msan  Energetics  for  a  Two-layer,  Free-surface,  Primitive- 
equation  Model  with  Open  Boundaries 

Consider  the  momentum  form  of  the  primitive  equations  for  a  two-layer  fluid 
with  a  free  surface: 


3  Wj 

at  +  'V7'Vkxf'V  -vpi 

+  (  ti  -  fi+1)/(phi)  +  AV2 \v. 
+  V  •  (hj  \Vj  +  h2  w2)  =  0 


3n2 

‘  (h2'v2>  =  0 


(Cl) 

(C2) 

(C3) 
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1 


where  i=l,2  and  and  n2  are  the  deviations  of  the  free  surface  and  the  inter 
face,  respectively  (see  Appendix  A).  Also, 

p1  =  gn1  9'  =  3(p2  "  pi^p 

P2  •  pi  -  g'foj  -  n2)  »i  =  Ti  1  +  Ti  j 

Note  in  (Cl)  that  we  have  used  the  traditional  form  of  Laplacian  friction. 

Now  define  the  kinetic  and  potential  energies  as 


=  bph,-  (u2  +  v?)  where  i =1 ,2 
P  =  Hp  (gn2  +  g'n2). 


(C4) 


Multiplying  (Cl)  by  ph^  (i=l,2)  and  using  (C2)  and  (C3)  we  obtain  the  kinetic 
energy  equations: 


aKl 

+  V  •  (  wiKi )  +  P9V  *  vVlnl^ 


2 

=  pgnjV  *  (fy  Wj)  +  pAhj  \vi  V  \vj 


and 


(C5) 


3K2 

3^-  +  ?  •  (av2K2)  +  pgv  •  (h2  w2nx) 

+  pg'v  ■  (h2w2n2)  =  pgnjV  •  (h2  w2)  (C6) 

+  pg'n2  v  •  (h2  w2)  +  pAh2  w2v2  v2 

Note  that  in  deriving  (C6)  we  have  assumed  g'Vrij  «  gVqj  and  f..=0  for  1*1, 2,3. 
The  potential  energy  equation  is  formed  by  multiplying  (C2)  by  pgnj  and  (C3) 
by  pg’n2  and  summing  the  results: 

+  pgnjV  •  (hj  Wj  +  h2  w2)  +  pg'n2v  •  (h2  w?)  =  0  (C7) 

Now  define  mean  and  perturbation  quantities  for  u,  v,  h,  and  n  such  that 

i  f  t  +  T 

(  )  =fj  0  (  )dt 


L 
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and 


(  )'  =  (  )  -  (  ). 

where  T  is  a  suitable  time  interval.  Also  define  the  kinetic  energy  of  the  mean 
flow  per  unit  area  as 


K.  =  >sphi  (u?  +  v/),  i = 1 ,2 

and  the  mean  kinetic  energy  of  the  eddy  flow  per  unit  area  as 
+  v^)  -  K. . 

Similarly,  for  the  potential  energy  of  the  mean  flow 

P  =  ^(gnj  +  g'n2  ) 

and  for  the  mean  potential  energy  of  the  eddy  flow 
P'  =  4p(gnJ^  +  g'?) 

The  equation  is  obtained  by  multiplying  (Cl)  with  i=l  by  ph^  w^: 
3K 

~  =  -  V  •  (  WjKj)  -  Yj  -  pgv  •  (n^  Wj) 

+  pgrijV  •  (h1  xv j)  +  pAhj  WjV2  \Vj 
where  for  layer  i 


(C8) 


(C9) 


(CIO) 


(Cll) 


(C12) 


Yi  =  *3p(ui2  +  v.2)V  •  (hi  wl )  +  ph^-tv  •  (u(\v!) 

-  ulV  •  \v[)  +  ph.v.(v  •  (vTvvT)  -  vjV  •  xvj). 
The  K[  equation  is  (C5)  -  (C12)  or 
9K* 

=  -  [V  •  \\T^  -  V  •  (  WjKj)]  +  Yj 

-  pgv  •  (njhj  v1  -  iijfij  xv^  +  pgir^v  •  (hj  xvj) 

-  njV  •  (\\l  xvj)]  +  pAlhj  WjT2  \v j  -  x^V2  wjl 
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Similarly,  for  K2  and  IC,  we  obtain 


3t~  =  "  v  *  (  vv2K2)  -  y2  -  pgv  •  (n^  \v2) 
+  pgr^V  •  (h2  vv2)  -  pg'v  •  (n2^2  \v2) 

+  pg'n2v  •  (h2  w2)  +  pAh,  ^v2  v7. 


2  w2v  w2 


=  -  [V  •  (  W2K2)  -  V  •  (  W21<2 ) ]  +  >2 

-  pgtv  •  (n1h2  w2)  -  v  •  (tT^FT^  \v2)i 


pg  tnjV  •  (h2  w2)  -  njV  •  (h2  v2)] 


-  pg'  [V  •  (n2h2  w2)  -  v  •  (npF,  \v,)j 


2"2  'V2I 


+  pg'  (n2v  (h2  w2)  -  n2v  •  {h2  w2)l 
+  pA[h2  -  h2  W2V2  vv2) 

The  P  equation  is  obtained  by  multiplying  (C2)  by  pgrij  and  (C3)  by 
and  summing  the  results:  1 

0P~  _  _ 

3t  =  -  pgnjV  .  (h:  Wj  +  h2  w2)  -  pgnjV  •  (hp]-  +  fi^) 

-  pg'n2v  •  (F2  vv2)  -  Pg'n2v  •  (hpvp 

By  subtracting  (C16)  from  (C7),  we  obtain  the  P'  equation: 

0p«  _ _ _ 

cTt  =  ‘  pg[n7v  ’  (hi  Vj  +  h2  \v2)  -  njV  •  (hj  vvj  +  h2  w2)] 

+  pgnjV  •  (h]  \vj  +  h^\v2) 

-  pg'  tn2v  •  (h2  w2)  -  n2v  •  (ff2  w2)j 
+  pg'n27  •  ( h2  w2) 


1 


When  the  terms  in  the  energy  equations  are  calculated  from  a  model  solution 
in  statistical  equilibrium  and  a  suitable  time  average  is  used,  the  tendency 
terms  are  negligible.  Regional  energy  budgets  can  be  obtained,  if  (C12)  -  (C17) 
are  also  Integrated  spatially.  Integrated  over  a  closed  basin  the  divergence 
terms  vanish,  but  in  an  open  domain  they  must  be  retained.  The  energy  balances 
are  conveniently  displayed  using  an  energy  box  diagram  such  as  Fig.  12a.  Each 
term  is  represented  by  an  arrow  in  or  out  of  the  box  for  the  energy  reservoir 
associated  with  that  equation.  When  identical  terms  with  opposite  sign  appear 
in  two  equations,  they  represent  a  conversion  of  energy  from  one  type  to  another. 
The  energy  transfers  shown  in  Fig.  12a  were  calculated  from  the  spatially  inte¬ 
grated  terms  in  C12)  -  (C17)  as  shown  below. 


Kj  terms 


{Kp  -  Kj} 


-  JJ V  •  (  WjKjJdx  dy 

{Kj  ->■  Kp  i^Yjdxdy 

{PW  Kj}  e  -  pgjjv  •  (njKj  \7pdx  dy 
{Tfj  P>  s  -  pg/JVp  •  (hj  \Vj)dx  dy 
{Kj  D}  e  \vp2  \Vj  dxdy 


Kj  terms 


{KjF  -  Kp  e  -JJ V  •  (  \vp^  -  vvppdxdy 
(Kj  -*■  Kp  e  JJy^dx  dy 

Kp  =  -  pg jj 7  •  (npj  Wj  -  npj 
•  -*  P' }  =  -  pg JJ mp  •  Ihj  \vp  -  np  •  (hx  \vpj 
’  -*■  0}  e  -  pljjlhl  \vp2w1  -  hj  \Vj  V2  Wjldxdy 


{PW 

{K 

{K 


)dx  dy 


dx  dy 


(C18) 


(C19) 
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Kg  terms 


{K2F  -►  k2j  =  #  •  (  w2K2)dx  dy 
{K2  -*•  KJj}  =  ^J*Y2dxdy 

{PW  -  k2j  =  -  #  •  [(gnj  +  g* n2)t»2  w2idxdy 
{K2  -*  F}  =  -  pjj (gnj  +  g 1  fT2 ) V  •  (h2  \v2)dxdy 
{K2  -*■  D}  =  -  p kjj h2  \v2V2  \v2dx  dy 


K2  terms 


{K£F  -*•  K2>  e  - JJv  •  (  w"2K2  -  \v2i<2)dxdy 

(K2  -*■  K£}  =  ^Y2dxdy  (C2‘ 

{PW  ■+  k2>  =  -  p  Jj\ 7  •  [(grrprTn^Th^w^  -  (gffj  +  g'n2)h2  \v2]dxdy 
{K2  -►  P* }  =  -  p JJ [(gnj  +  g'n2)v  •  (h2  \v2)  -  (gr^  +  g'n2)V  •  (h2  \v2)]dxdy 
{Kg  -*■  0}  =  -  pA  jj [h2  ^V2  w2  ■  ^2  ^2^  v^21  dx  dy 

F  terms 

{Kj  -*■  P)  =  -  pg JJ HjV  •  (hj  \Vj)dxdy 

{F  -*■  K2>  =  pjj (gr^  +  g'n2)v  •  (h2  \v2)dx  dy  (C2; 

{F  -*•  P' }  =  pj^fgn^  •  (hj'w^  +  h2  w2)  +  g'n2V  •  (h^w^)ldxdy 


P1  terms 


{Kj  ->  P' }  =  -  pg//  hjV  •  (hj  Wj)  -  njV  •  (hj  Wj)]  dx  dy 

//  {C: 
{P*  ■>  K£}  s  pj^gni  +  g'n2;v  •  (h2w2)  -  (gnj  +  g'n2)V  *  (h2w2))dxdy 

{F  -»•  P'}  =  pff  [g^V  •  (h[\v[  +  h2  w£)  +  g'n2^  •  (h2  wpidx  dy 


The  bracket  notation  is  dropped  in  Section  8.  In  our  results  |Kj  -»  4 1  » 

I4F  -*■  K|  +  PW  ■*  Kj|.  Thus,  any  controversy  over  the  formulation  of  these  terms 
(Harrison  and  Robinson,  1978)  should  not  cloud  the  interpretation  of  the  basic 
results  in  Fig.  12.  Also,  the  difference  between  the  frictional  formulation  in 
the  models  and  the  energetics  did  not  result  in  any  serious  imbalances  in  the 
energy  equations. 
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